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We report detailed systematic measurements of the spatial variation in electronic states in the high 
Tc superconductor La2-xSrj;Cu04 (0.04 < x < 0.16) using "^^Cu NQR for ''^Cu isotope enriched 
poly-crystalline samples. We demonstrate that the spatial variation in local hole concentration 
^'^ xiocaii"^ x) given by ^'^xiocai = x Axiocai, where x is the nominal hole concentration and 
^^Axiocai is defined as the amplitude (or extent) of the spatial variation, is refiected in the frequency 
dependence of the spin-lattice relaxation rate ®^1/Ti across the inhomogeneous linebroadening of 
the "^^Cu NQR spectrum [P.M. Singer et al, Phys. Rev. Lett. 88, 47602 (2002)]. We show that 
compared to nominal x, the electronic state in certain regions of the Cu02 plane are locally more 
metallic {^'^xiocai = x +^''^ Axiocai) while others are more insulating i^^xiocai = x — ®^ Axiocai)- By 
using high precision measurements of the temperature dependence of ^^l/Ti at various positions 
across the "^^Cu NQR lineshape, we demonstrate that Axiocaii"^ 0) increases below 500 - 600 K 
and reaches values as large as ^^Axiocai/x ~ 0.5 in the temperature region > 150 K. We find a 
substantial overlap between the ®^Cu NQR spectrum of samples with different x, and find that 
the extent of the overlap increases with decreasing temperature. By incorporating the random 
positioning of """^Sr donor ions in the lattice in a novel approach, a lower bound to the length scale 
of the spatial variation ^^Rpatch is deduced by fitting the entire ®^Cu NQR spectrum (including the 
"B" -line originating from ®"^Cu sites with "'"^Sr ions directly above or below) using a patch- by-patch 
distribution of the spatial variation ^'^xiocai with the patch radius ^^Rpatch 3.0 nm (= 8-10 a, 
where a is the lattice spacing) as the only free parameter. A corresponding upper bound to the 
amplitude of the spatial variation Axpatc/i(oc 'i/^^Rpatch) is deduced within the patch- by-patch 
model, and consistent results are found with ^'^ Axiocai determined from the frequency dependence 
in ^^l/Ti. Using our pool of ®^Cu NQR data, we also deduce the onset temperature Tq{>^ 400 K) 
of local orthorhombic lattice distortions which, in the region x > 0.04, is found to be larger than 
the onset temperature of long range structural order. 



PACS numbers; 74.80.-g, 74.72.-h, 76.60.-k 
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I. INTRODUCTION 

NQR (nuclear quadrupole resonance) and NMR (nu- 
clear magnetic resonance) have proven to give unique in- 
formation regarding the inhomogeneous electronic state 
in the Cu02 plane of various high cuprates 1, 2, \ ^ 
SmSSIlfllllllllll including La2-x[Sr,Ba],Cu04, 
[La,Nd,Eu]2-a;[Sr,Ba]a;Cu04 and La2Cu04+5. In these 
materials, the high Tc superconductivity is achieved 
by doping holes into the Cu02 plane. In the case 
of La2-a;[Sr,Ba]a:Cu04 and its Nd or Eu co-doped 
compound [La,Nd,Eu]2-a;[Sr,Ba]a;Cu04, hole doping is 
achieved by substituting ions of different ionicity, thereby 
creating an alloy with intrinsic inhomogeneities caused by 
chemical substitution Over the past several years, 

studies using ^^Cu NQR and NMR wipeout 0,11 HlH 
ITU have characterised the glassy nature of the slowing 
down of the stripe inhomogeneity in these materi- 
als, where the Coulomb potential from the distorted lat- 



tice slows down spin and charge density waves. In the 
case of La2Cu04+5, hole doping is achieved by super- 
oxygenation, and it has been shown that the high mo- 
bility of the excess oxygen atoms results in electronic 
phase separation [ la. IitI lla | between the superconduct- 
ing and antiferromagnetic phase, as evidenced by ^'^^La 
NMR nil]. 

No clear picture has emerged which discerns and 
relates the effects of genuine electronic phase separa- 
tion, stripe modulation, and the random substitution 
of donor ions. On the other hand, recent STM (scan- 
ning tunnelling microscopy) studies on the surface state 
of Bi2Sr2CaCu208+5 cleaved at low temperature 19] re- 
veal spatial variations of the electronic state on a short 
length scale ~nm in the surface plane. Whether such 
nm modulations are universally observable in the bulk 
and other high Tc cuprates remains to be seen, but the 
STM results have enhanced the interest and potential 
impact of the spatial inhomogeniety of the electronic 
properties in cuprates. One emerging counter example 
is YBa2Cu30j, where NMR [2y] measurements in- 
dicate that the spatial inhomogeneity in YBa2Cu30j, is 
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less than in Bi2Sr2CaCu208+5 or La2_2;Sr2;Cu04. 

Unlike scattering techniques which probe coherent phe- 
nomena over length scales larger than tens of nm's, NMR 
and NQR are strictly local probes which make them an 
ideal tool for revealing the short length scale physics of 
the Cu02 plane. Among earlier reports of the inhomo- 
geneous state at short length scales in La2-KSr^Cu04 
are: the sphtting of the "^^Cu NQR and NMR lines due 
to inequivalent Cu sites known as the 'A" and "B" -sites 
resulting from different EFG (electric field gradient) ten- 
sors 1^, JJ due to the presence of nearby '*'^Sr ions, the 
drastic broadening of the zero field ^^Cu line at tempera- 
tures below ~ 4 K Hill, the '^^cu NMR line broadening 
from short length scale modulations in orbital shifts 0, 
and also the ^^Cu NQR and NMR wipeout 0,11 [13 as 
a result of the glassy slowing down of the stripe inhomo- 
geneity in the temperature region < 100 K. 

These evidences clearly support the existence of some 
sort of short length scale inhomogeneity in the Cu02 
planes of La2_a;Sra;Cu04 which raises questions regard- 
ing recent theoretical debates of a "universal electronic 
phase diagram", including La2-2;Sra;Cu04 1213 . which 
are based on the assumption that hole doping is ho- 
mogeneous. More recently, we reported evidence of an 
inhomogeneous electronic state in La2-2:Sr2:Cu04 using 
^^Cu NQR in ^'^Cu isotope enriched samples in the range 
0.04 < a; < 0.16 [13. As discussed in Ref. [12, the 
two essential ingredients for the spatial variation in local 
hole concentration ^^xiocaiiy^ 0) are (a) the fact that hole 
doping is achieved by random substitution of donor ions 
in the lattice with different ionicity, and (b) the presence 
of a short electronic length scale ^^Rpatch ^ 3.0 — 4.0 nm 
for the spatial variation. To the best of our knowledge, 
our results reported in this paper and in Ref. jl2j are the 
first of its kind to detect the temperature dependence of 
the inhomogeneous electronic state in La2_a;Sra;Cu04, or 
any other high Tc materials with quenched disorder. 

Most recently, ^^O NMR in high-quality 
La2-a;Sr2;Cu04 single crystals in the range 
0.035 < X < 0.15 also reveal a substantial inhomo- 
geneity in the local electronic state The spatial 
variation in local hole concentration ^"^xiocai determined 
using ^^O NMR is found to reveal a consistent value 
with ^^xiocai reported here and in Ref. |0. The spatial 
variation ^"^xiocai E3 is determined through the observed 
spatial variation in the spin susceptibility across the 
Cu02 plane, whereas ^^xiocai is determined through the 
spatial variation in the EFG which is purely a charge 
effect. The fact that ^^xiocai — ^"^ xiocai therefore implies 
that the inhomogeneity in the spin and charge channels 
are highly correlated. 

As presented in Ref. Ts'l, we also compare the ex- 
tent of the inhomogeneity in high-quality single crystals 
[T3I I and our poly-crystalline samples (reported here and 
in Ref. ^3) comparing the extent of the frequency 
dependence of across the ^"^Cu NQR spectrum. 

We find consistent results between single crystal and 
poly-crystalline samples [l3j |. which immediately estab- 



lishes that the inhomogeneity in the electronic state is 
an intrinsic phenomenon in La2-a;Sra;Cu04. It also rules 
out certain claims that our poly-crystalline samples are 
somehow more inhomogeneous than those reported else- 
where. Such claims on the quality of our poly-crystalline 
samples originated as an excuse to negate the fact that 
hole doping is inhomogeneous in La2-a:Sra:Cu04 T^. In 
this paper we report the details of the ®'^Cu NQR results 
and analysis originally outlined in Ref. |l2l. and while 
the ^^O NMR results are left to Ref. [13, we mention 
the conclusions of the "'^'^O NMR results wherever appro- 
priate. 

The dynamics of the electronic state are probed using 
the spin-lattice relaxation rate is a mea- 

sure of the local spin fluctuations in the Cu02 plane at 
the ^^Cu NQR frequency ^^vq (ranging between 32 — 40 
MHz). Due to the experimental ease of the measurement, 
the temperature dependence of ^^1/Ti is generally only 
measured at the peak of the spectrum. The most striking 
feature of our ^^Cu NQR data, however, is that ^^1/Ti 
shows qualitatively different temperature dependence de- 
pending on where along the NQR line it is measured [^ . 
This implies, without using any kind of model, that cer- 
tain regions of the Cu02 plane are more metallic while 
others are more insulating • 

In order to measure the frequency dependence in 
^^1/Ti and obtain any quantitative information over a 
wide temperature range, it is essential to separate the 
signal from ^^Cu and ^^Cu isotopes. Previous work by 
S. Fujiyama et al. [H measured (63,65) ^^^^ for natu- 
rally abundant Cu in which they successfully deduced 
that a substantial frequency dependence exists across the 
(63,65) NQj^ spectrum, however, no statement beyond 
that could be inferred. Through ^^Cu isotope enrich- 
ment, together with a systematic study as a function of 
X, we can use the extent of the frequency dependence in 
^^1/Ti across the ^^Cu NQR line to determine the local 
hole concentration ^^xiocai given by 

'^^Xlocal^X±^^ /^XlocaU (1) 

where ^^Axioca; is the characteristic amplitude or extent 
of the spatial variation ^^xiocai- 

In Fig. ^ we show an example of the frequency de- 
pendence in 63i/r^ across the ^^Cn NQR Hne at 300 K. 
We first remark on the overlap between samples with 
different nominal hole concentration x, shown in Fig. 
n^b). Taking x = 0.115 as an example, we see that the 
upper (lower) half intensity point of the x = 0.115 A- 
line (or B-line) roughly coincides with the peaks of the 
X — 0.16(0.07), which implies that the characteristic lo- 
cal hole concentration xiocai is given by 

•^local 

< 0.16 

in the more metallic regions, and by ^^xiocai > 0.07 
in the more insulating regions. Using these limits and 
Eq. QJ, we estimate the upper bound to the amplitude 
^^Axiocai < 0.045 for x = 0.115 at 300 K. The reason 
that only an upper boundary to ^^xiocai (and an upper 
bound to ^^Axiocai) can be deduced is that we have as- 
sumed that the linewidths in Fig. are dominated by 
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^^xiocai{^ x) alone. In fact, as we calculate in section IV, 
there is a substantial intrinsic lattice linebroadening (de- 
fined as /^viatt) which is independent of the broadening 
arising from ^^xiocai- ^i^iatt originates from the distribu- 
tion in local EFG values due to the random positioning 
of +^Sr ions in the lattice, and we calculate its size by 
using a point charge lattice summation. In Fig. |21 we 
illustrate the extent of Aviatt as the dashed curve, which 
is comparable to the experimentally observed linebroad- 
ening. In order to determine a better estimate of ^'^xiocai 
in a model independent way, we must therefore measure 
^^l/Ti at various positions across the line as shown in 
Fig. IJa) and Fig. 

In order to determine ^^xiocai from we first de- 

fine the values ^^1/Ti,a and ^^1/Ti,b taken at various 
positions across the NQR lineshape (shown in Fig. 
63l/r{°] defines measured at the CG (center of 

gravity) of the A-line, defines ^'^l/Ti measured 

at the half intensity of the upper (-I-) frequency side of 
the A-line, 2 ^'^ the lower half intensity side, and 

^^l/T^^^^^"-' at the lower one-tenth intensity of the A- 
line. We also define the corresponding quantities for the 
B-line. We shall use these definitions throughout this pa- 
per to deduce ^^xiocai as a function of temperature and 
X. Next, we assume that *^^l/r|°j and ^^l/T^% taken 
at the CG of each spectrum represents the value 
for the nominal hole concentration x, regardless of the 
underlying linebroadening. We can then compare 
at the upper and lower half intensity points of the A-line 
(or B-line) to nominal values ^^1/T^^\ from different x. 
As we shall show in detail in section III, ^^l/T^^J for 
X = 0.115 at 300 K is the same as ^^l/T^^l for x = 0.08, 

while ^^l/T^^2 for ^ = 0-115 is the same as ^^l/r/°i 
for X = 0.15. Using Eq. ^ we therefore deduce that 
^^Axiocai = 0.035 for x = 0.115 at 300 K. Since we do 
not have the continuous function as a function 

of nominal x (i.e. we have a limited number of sam- 
ples), we interpolate ^^1/t|°4 between the discrete set of 
X we do have. The same argument applies to the B-line, 
and we shall demonstrate that ^^Axiocai deduced using 
the B-line is consistent with that from the A-line. In 
order to be consistent with Ref . JJJ , we define the char- 
acteristic amplitude ^^Axiocai with respect to the half 
intensity points, however, following the same proceedure 
using 1/T^~^^^^'' instead results in an overall 20 - 40 % 
increase in ^^Axiocai and does not effect our conclusions. 
Using the frequency dependence of ^"^l/Ti, we therefore 
deduce that within each sample with nominal x there ex- 
ist certain regions in the Cu02 plane with higher and 
lower local hole concentrations {^^xiocai 7^ x). System- 
atic studies as a function of temperature then allow us 
to follow the temperature dependence of ^^xiocai- 

The next task in this paper is the model analysis of 
the NQR spectrum. For this task, we develop a static 



real space picture consisting of a patch-by-patch spatial 
variation for ^^xiocai [^]- Using a point charge lattice 
summation of the EFG together with a patch-by-patch 
spatial variation, we successfully fit the NQR spectrum 
(including the spectral overlap for different x and the B- 
line) with one free parameter {^^Rpatch)- As detailed in 
section IV, the simulation of the EFG consists of two es- 
sential ingredients: (a) the random positioning of '''^Sr 
donor ions in the lattice, and (b) the presence of a short 
length scale ^^Rpatch ^ 3 nm which describes the in- 
homogeneous electronic state. In our model, the Cu02 
plane is sectioned into patches (i.e. circles) of equal ra- 
dius Rpatchi— 3 nm), and each patch i has a different 
local hole concentration x]^^^i{^ x). The value of x]^^^i 
is determined by the random number of +^Sr donor ions 
in the vicinity of patch i. The holes are placed uniformily 
within patch i, and the resonance frequency of the ®^Cu 
nuclei within patch i is calculated in a self-consistent way 
using our pool of NQR data. By sampling a large num- 
ber of patches (~ lO'') and calculating the resonance fre- 
quency for each patch, we build up a histogram spectrum 
of resonance frequencies which constitutes the simulated 
NQR lineshape. ^^Rpatch is the only free parameter in 
the simulation, and we optimize ^^Rpatch to fit the NQR 
spectrum, as shown in Fig. [3 

Using ^^Rpatch = c» is equivalent to uniform homoge- 
neous doping where ^^xiocai = x throughout the plane 
(i.e. ^^Axiocai = 0). This is also equivalent to excluding 
ingredient (b) listed above and calculating the broaden- 
ing solely from the random positioning of +^Sr donor 
ions in the lattice (i.e. Ai^iatt)- Clearly, however, in the 
case ^^Rpatch = 00 the observed linewidth is underes- 
timated by Aiyiatt alone, likewise for the B-line. Using 
a very short length scale ^^Rpatch = 1-6 nm overesti- 
mates the NQR linebroadening, whereas the optimum 
value ^^Rpatch — 2.6 nm results in the best fit to the 
data, including the B-line. In Ref. T^l we also show 
that ^^Rpatch thus deduced is consistent with ^"^ Rpatch 
determined in a similar way using the inhomogeneous 
linebroadening of the ^^O NMR spectra. 

Once ^^Rpatch is optimized to fit the NQR spectra, we 
deduce the spatial variation ^^Xpatch given by 

Xpatch — X ziz AXpdicji^ (2) 

where Axpatch is the calculated amplitude for ^^Xpatch- 
^^Axpatch is found to follow the simple relation 
^^Axpatch oc 1/^^Rpatch- The results of the EFG simu- 
lation {^^Xpatch) are found to be consistent with ^^xiocai 
deduced from the frequency dependence in Since 
^^xiocai is determined in a model independent way, the 
fact that ^^Xpatch —^^xiocai justifies the use of our patch- 
by-patch model. We note, however, that apart from the 
linebroadening mechanisms (a) and (b) discussed above, 
we ignore all other possible NQR linebroadening mech- 
anisms in our model. Other possible mechanisms in- 
clude distributions in local lattice distortions (section 
V). Our calculation in section IV ignores such possi- 
bilities which implies that ^^Rpatch is a lower bound. 
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that ^^Axpatch{oc l/^^Rpatch) is correspondingly an up- 
per bound, and therefore that ^^Xpatch is an upper bound- 
ary. The fact that ^^Xpatch ^ ^^xiocai, however, indicates 
that our simulation of the NQR spectra and the patch- 
by-patch model is appropriate. 

In section V we introduce the effects of the local or- 
thorhombic distortions into the simulation. We find that 
the NQR data separates into two temperature regions, 
above and below Tq{> 400) K. For T > Tq, our data in- 
dicates that the lattice is locally in the tetragonal phase 
and we use this region to estimate ^^Rpatch and ^^Xpatch 
as discussed above. For T < Tq, however, we find ev- 
idence for local orthorhombic distortions which, in the 
case oi X > 0.04, start somewhat above the structural 
phase transition temperature Tst(< 515 K) according to 
neutron diffraction studies ^22] . In the doping region 
X > 0.04, we show the systematic differences between 
the local and long range orthorhombic order. 

The rest of this paper is sectioned as follows: in section 
II we describe the experimental characteristics including 
the sample characteristics and NQR background. In sec- 
tion III we present the experimental results which include 
the frequency dependence of and the determina- 

tion of ^^xiocai- We also present the temperature depen- 
dence of the NQR spectra. In section IV we describe 
the analysis of the NQR spectra using a patch-by-patch 
model for the spatial variation in ^^xiocai- We present 
the details of the EFG simulation incorporating random- 
ness effects, and we optimize the length scale ^^Rpatch 
(and ^^Xpatch) to fit the experimental data. In section 
V we incorporate the local orthorhombic distortions into 
the EFG simulation. In section VI we then present our 
conclusions. 



II. EXPERIMENTAL CHARACTERISTICS 

A. Samples 

All of the single phase, poly-crystalline samples of 
La2-2;Sr3;Cu04 used in this study were prepared using 
conventional solid state reactions 23s3 '2M ■ We mix pre- 
dried LasOg (99.995%), SrCOg (99.995%), and ^^CuO 
(99.995%) with correct nominal compositions using an 
agate mortar and pestle until an intimate mixture is ob- 
tained. A pre-reaction is carried out for 20 h in a box 
furnace at 850°C followed by repeated grindings and sin- 
terings (also 20 h) at temperatures between 950° C and 
1000°C. During these initial firings the materials are 
made into low density pellets with a hand press. The 
samples are then pelletized with the pneumatic press and 
high temperature annealed in flowing O2 gas at 1100°C 
to 1150°C for 24 h to 48 h before a slow and controlled 
cooling cycle that includes low temperature annealing at 
800°C (24 h) and 500°C (24 h). The long annealing in 
O2 insures that the oxygen content is uniform and sto- 
ichiometric and high annealing temperatures insure fast 
reaction kinetics. By using a large number of grindings 



(typically 5-8), we achieve high quality poly-crystalline 
samples. 

The most important test for the quality of our samples 
was discussed earlier (and presented in Ref. |0) where 
we confirmed that the frequency dependence of 
across the ^^Cu NQR spectrum in our poly-crystaUine 
samples is the same as in high-quality single crystals 
grown m Tokyo [3. This confirms that the inhomo- 
geneity we report in this paper and in Ref.'s are 
intrinsic to La2_a;Sra;Cu04. Apart from ^^Cu NQR itself, 
we characterise our poly-crystalline samples by measur- 
ing the superconducting transition and the room tem- 
perature lattice parameters. We deduce by measuring 
the Meissner signal using a SQUID magnetometer in the 
field cooled mode with a constant applied magnetic field 
of 15 Oe , the results of which are shown in Fig. |31 We 
find that the onset temperature of the diamagnetic sus- 
ceptibility, Tc, and the volume fraction agree with previ- 
ous measurments in poly -cr y stalline samples prepared in 
a similar fashion [23 . |24L l25l | . We note that the onset of 
the diamagnetic susceptibility is somewhat tailed in the 
case of X = 0.20, consistent with 2i,2A,2M- In Ref. 
it is even argued that bulk superconductivity disappears 
in the overdoped region, which coincides with the disap- 
pearance of LRO (long range order) into the orthorhom- 
bic structural phase. As we present in section III, we 
also find a qualitative change in in the overdoped 

region. For these reasons we limit our determination of 
^^xiocai in section III to x < 0.16. 

We deduce the lattice parameters with an X-ray 
diffractometer using the Cu Ka^ line. The lattice pa- 
rameters presented in Fig. 0] (the values of which we use 
in the analysis in section IV) show close agreement with 
previous results [^. In this paper, we define and bo 
to be the in-plane lattice parameters of the orthorhombic 
cell and [ao,bo,c] as the vector along the orthorhombic 
directions. Likewise, we define a and b{= a) as the lat- 
tice parameters in the tetragonal cell and [a, 6, c] as the 
vector along the tetragonal direction. Our typical X-ray 
linewidths were comparable to the instrumental resolu- 
tion with a HWHMoi Sinstr ~ 0.025° determined using a 
high quality Si standard. Previous higher resolution X- 
ray diffraction experiments '25^ (where Sinstr ~ 0.01°) 
were able to use the observed X-ray linewidths and 
overlaps between different samples x to deduce a lower 
bound to the compositional distribution with a HWHM 
of 5xsr ~ 0.01. Only a lower bound to 5xsr could 
be inferred because X-ray diffraction takes a spatial av- 
erage over length scales larger than tens of nm's '2^. 
They attributed 5xsr ~ 0.01 to imperfect mixing inher- 
ent to the solid state reaction. Similar conclusions were 
reached in Ref. 0| using the overlap of X-ray diffrac- 
tion and neutron Bragg peaks between different x, where 
5xsr ~ 0.015 was deduced. 

We confirmed the same order lower bound 5xsr — 
0.013 in single crystal a; = 0.15 using EMPA (electron 
micro-probe analysis) whose length scale is determined 
by the focus area of the beam ~ 1 /im. This means that 
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even high-quality single crystals whose spatial variation 
of +^Sr content is as small as Sxsr — 0.013 {<Si^^ ^xiocai) 
over 1 /im, exhibit a large variation in the local elec- 
tronic states. Furthermore, we will demonstrate in sec- 
tion III that A.xiocai is temperature dependent and in- 
creases with decreasing temperature. This means that 
the variation in +^Sr content from the solid state reac- 
tion alone cannot account for our findings. Instead, a 
thermodynamic process must be involved. 



B. NQR Background 

The essential condition for NQR is that the ground 
stat e sp in "/ of the nuclear isotope n satisfy "/ > 1 
[2^ l27j . In the case of Cu, there exist two stable iso- 
topes n = 63(65) with natural abundances of 69 (31) % 
respsectively, both with spin (^'^-^s)/ — 3/2. This is both 
an advantage and a disadvantage. It is generally a dis- 
advantage in the sense that their quadrupole moments 
"■Q happen to be similar ^^Q/^^Q = 1.081, and therefore 
their resonance spectrum tend to coincide. This is the 
case in La2_a;Sra;Cu04 for x > 0.05 where the NQR spec- 
trum from each isotope is broad enough that they merge 
in the temperature range of interest, and therefore ^^Cu 
isotope enrichment is essential to separate multiple lines. 

The observed NQR frequency ^^i^q is given by the total 

EFG (electric field gradient) tensor v/j^''^'' (where P and 
7 are any two orthogonal spatial directions) surrounding 
the resonating ^^Cu nucleus, and its quadrupole coupling 
constant ^^Q. In terms of magnetic eigenstates m/, NQR 
for 1 — 3/2 corresponds to transitions between mj = 
I ± 1/2| and = | ± 3/2| eigenstates 26, 27]. The EFG 
cannot split the ± degeneracy of the magnetic eigenstates 
since it is of charge origin. The most general expression 
for the NQR tensor is given bv [2ll27|: 



63, (/3,7) 
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T/(/3,7) 
''tot I 



(3) 



where v/j^'''^'' is the traceless symmetric EFG tensor. 
Vfgf'^^ can be rotated and diagonalised to point towards 
the principle axes a of the EFG. This leads to the fol- 
lowing relation 
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In the tetragonal phase, the principle axes a are along 
the crystal axes [a, 6, c], and satisfy the general traceless 
condition V'^ + + V°' = 0. The direction of the largest 
component IV] (where IV] > \V^\ > \V°-\) defines the 
main principle axis in the tetragonal phase. The asym- 
metry parameter rj (0 < 1] < 1), defined as 



ya _ yb 



(5) 



is typical ly s mall 77 < 0.06 for the planar Cu site in the 
cuprates and for both A and B-sites [l^. The ob- 
served NQR frequency ^^vq is then given by 
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•F,:,.(l + W3) 



1/2 



(6) 



In the present case, r] is small [r] < 0.06) and to a very 
good approximation 
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The value of ^^Q used varies between groups, but is 
generally taken to be either ^^Q — —0.211 barns |30| 
which is based on theoretical calculation, or ^^Q = —0.16 
barns jU which is based on ESR measurements in 
cuprate salts. Since we are dealing with experimental 
results, we choose to use ^^Q — —0.16 barns [sj. As 
discussed in appendix A, however, our final results of the 
length scale ^^Rpatch is insensitive to the absolute value 
of '^^Q. 

In NQR, the direction along which the spin-lattice re- 
laxation rate ^^1/Ti is measured is given by the main 
principle axis of the EFG (i.e. the c-axis), and cannot 
be changed externally. The most general expression for 
^^1/Ti in NQR with relaxation by a single magnetic pro- 
cess is given by [s^ 
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63„,2 



|2 x'l (q,w„) 



(8) 



.63 32 _ 40 MHz) is the NQR fre- 

11.285 MHz/Tesla is the nuclear gyro- 



where a;„/27r 
quency, ^^^/2-k 
magnetic ratio, and q is the reciprocal lattice vector (in 
tetragonal notation) . 

^^A{(\)-^ is the hyperfinc form factor perpendicular to 
the main principle axis. The two perpendicular directions 
(_L) in the present case both lie in the Cu02 plane (i.e. 
in the a — h plane). The q dependence of the in-plane 
value 63^(q)-L ig given by 

^^A(q)^ = A^ + 2B- {coa{q^a) + coa{qya)) , (9) 

where A± ~ 38 kOe/ iib and i3 ~ 42 kOe/jiB are 
found to be the same in La2-sSr2,Cu04 [s^, l34| and 
YBa2Cu30y [33 compounds, x'i(qj'^ra) (in emu/[mol 
f.u.]) is the in-plane imaginary part of the dynamic 
spin susceptibilty of the electrons. According to Eq. 
0, there is a large contribution to from the 

anit-ferromagnetic wavevector Oaf = (7r/a,7r/a) where 
63A(Qaf)^ = -139 kOe/nB iHl. Since xT(q,t^n) is 
peaked in the vicinity of q = Q af , this implies that 
^^1/Ti is a local probe of the low-frequency (i.e. hujn <C 
ksT) in-plane anti-ferromagnetic spin fluctuations. ^^7^ 
enters into Eq. © because the spin-lattice relaxation 
process is magnetic in origin, whereas the splitting of the 
energy levels ujn is determined by the EFG which is of 
charge origin. 
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We measure by applying the following pulse 

sequence 



-t 



echo. 



(10) 



where an r.f. pulse tt inverts the nuclear magnetization 
M{t) at time t prior to the spin-echo sequence in brackets. 
The time dependence of M{t) is determined by taking the 
integral of the echo in the time domain as a function of 
delay time t. ^^\/Ti is then determined by fitting M{t) 
to the standard recovery form 



The " static" time-scale is given by the spin-echo refocus- 
ing time 2r — 24 /is (Eq. l(Tn)l'). 2t defines the effec- 
tive exposure time (or equivalently, an effective shutter 
speed of 1/(2t)) for the NQR experiment. All dynamic 
effects with time-scales shorter than 2t are unobservable 
through the NQR spectrum. 



III. EXPERIMENTAL RESULTS 



'l/Ti Results 



M{t) = M{oo) + {M{Q) - M{oo))- 



exp 



Tl 



, (11) 



Before presenting the frequency dependence of 
we first present the temperature dependence of ^'^1/Tf\ 



appropriate for NQR with "/ = 3/2. 

The results of ^^l/Ti were independent of the r.f. pulse 
width TT in Eq. In the case of a: < 0.16, M{t) fit weU 

to Eq. Hll() , implying that the spin-lattice relaxation rate 
is dominated by a single relaxation mechanism. In the 
case of X > 0.16 below < 100 K, however, detailed mea- 
surements of M{t) revealed a multi-exponential recovery, 
implying that ^•^l/Ti is distributed at a fixed frequency. 
In appendix B we use the multi-exponential recovery of 
M{t) to deduce the intrinsic lattice broadening ^vf^^^ of 
the NQR line. In the case of a; = 0.16, we find that 
^^Tatt ~ 0.62(±0.07) MHz, which is consistent with our 
independent results of the EFG simulation Ai^jatt = 0.49 
MHz presented in section IV. All ^^\/Ti data presented 
in this paper (except for appendix B) correspond to re- 
sults of force-fits (in the case x > 0.16) deduced using Eq. 
()ll|l . The force-fit values from Eq. Ullfl correspond to 
the average value of the underlying distribution in spin- 
lattice relaxation rate. 

As a result of the distribution in the spin-lattice re- 
laxation rate, the value of ^■^1/Ti determined using Eq. 
((TT|l showed a slight dependence on the pulse separation 
time T in Eq. H10|l . The most likely explanation is that 
using a long r results in suppressed signal intensity from 
^^Cu nuclei with fast spin-spin relaxation rates {^^1/T2). 
Since ^'^Cu nuclei with fast also have fast compo- 

nent in ^^1/T2 [33, one would expect a decrease in the 
average value of ^^1/Ti with increasing t. In the case of 
X = 0.16, increasing r from 12 to 24 /is resulted in a uni- 
form ~ 10 % decrease for ^^l/Ti across the NQR line. In 
order to retain the maximum possible signal intensity and 
measure the most representative value of ^^1/Ti at each 
frequency, we used the shortest possible value r = 12 ^s 
the experiment allows for all measurements in this paper. 
Using the same experimental conditions throughout also 
allowed for a consistent determination of ^'^xiocai- 

Finally, we note that there are two relevant time-scales 
in the NQR experiment. The first is the "dynamic" time- 
scale associated with Since ^'^1/Ti measures the 
spin fluctuations specifically at the NQR frequency ^^vq 
(Eq. ©), the appropriate "dynamic" time-scale is given 
by ~ 0.027 /iS. The other relevant time-scale 
is the "static" one associated with the NQR spectrum. 



and ^"^I/tI"^ taken at the CG of each spectrum. In 
Fig. 13 and Fig. |^ we present the temperature depen- 
dence of ^^1/t[°)^ and ^^1/t[^1, respectively. While the 

solid curves through the ^^1/Tf\ data in Fig. |S1 are 
guides for the eye, the solid curves in Fig. |B1 represent 

U-^'^1/t[°}^ data taken from Fig. [3 We plot Fig. 
|6| as such to directly compare the temperature depen- 
dence of -^^I/tI"]) with ^H/t[% where e is a uni- 
form scaling factor taken to be e = [0.87,0.87,0.90,0.84] 
for x = [0.20, 0.16, 0.115, 0.07], respectively. As shown in 
Fig. 1^ we find semi-quantitatively the same temperature 
dependence between the CG of the A and B-lines, with 
an overall 10 - 16 % smaller value of ®^1/Ti.s compared 
to ^3l/Ti,A, i-e. ^^Ti^a/^'^Ti^b = 0.9-0.84. Previous re- 
ports in non-isotope enriched samples 38, 39] also found 
that ^^Ti^a/^^Ti^b < 1- 

At first glance of Fig. El our data for x = 0.16 and 
X ~ 0.20 below 150 K may suggest that the ratio 
^^Ti^a/^^Ti,b tends to decrease - 10 - 15 % with de- 
creasing temperature, similar to previous trends reported 
in 38] for x = 0.20. We point out, however, that cau- 
tion must be used when deducing any significance from 
temperature variations in ^^Ti a/^^Ti b- The reason for 
this stems from the fact that spin-lattice relaxation rate 
is distributed at a fixed frequency, and that finite r = 12 
/iS conditions are used. As discussed in section II, the 
force-fit value will tend to underestimate the CG 

of the underlying distribution, the extent of which will 
depend on t and the one-to-one correspondence between 
the spin-spin relaxation rate and spin-lattice relaxation 
rate for all ^"^Cu nuclear sites. However, it is known that 
the spin-spin relaxation rate shows a complicated tem- 
perature dependence for both A and B-sites. In partic- 
ular, the Gaussian-like component of the spin-spin re- 
laxation rate (^'^ 1/720 40]) tends to get smaller with 
decreasing temperature while the Lorenztian- like com- 
ponent (^^1/T'2l) becomes dominant. Below the onset 
of ^^Cu wipeout T]\jQji, the spin-spin relaxation rate be- 
comes totally Lorentzian-like [a, S B uM both A and 
B-sites. Slight differences in the temperature dependence 
of *^^1/T2G and '^^I/Tsl between the A and B-sites com- 
plicate the one-to-one correspondence with spin-lattice 
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relaxation rate, and result in different underestimations 
of between the A and B-sites with decreasing tem- 

perature. Such factors must be considered when inter- 
preting the temperature dependence in ^^Ti a/^^Ti^b for 
< 150 K and a; > 0.16. 

In Fig. [7| we present the temperature dependence of 
at various frequencies across the line, labeled 
according to Fig. El In the background of Fig. [3 
we also show ^^l/T^^^ taken from Fig. |S1 which repre- 
sent curves of constant x. The most surprising discov- 



ery of the present work is that ^^1/T^~^\ ' and 

^^1/t| ^/^"^ all show qualitatively dijferent teurperature 



.(-) 



dependence. For example, ^'^l/T^^^ for x — 0.07 exhibits 
semi-quantitatively the same behaviour as for 
X = 0.115, while ^^1/T^''> for x = 0.07 exhibits semi- 



quantitatively the same behaviour as ^^l/T^]^ 



for X 



0.04. This is consistent with the fact that the up- 
per (lower) frequency side of the NQR spectrum for 
X = 0.07 roughly coincides with the peak NQR frequency 
of a: = 0.115 {x ~ 0.04) respectively, as shown in Fig. Q 
We also note that with decreasing temperature, 



63 1 



.(-1/10) 



.(0) 



do not exactly follow ^n/T|' 



for a given x, which indicates that °'^Axiocai is gradually 
growing with decreasing temperature. 

In Fig. |Slwe illustrate the details of the process used 
to extract ^^xiocai from the data in Fig. [7| We first take 
the data from Fig. □ and plot 



(+) 63 
A ' 



and 



^'^l/T^^\ for each a; at a fixed temperature of 300 K. For 
clarity we connect the ^^1/Ti,a data for fixed x by the 
dashed black lines. We then create a smooth interpola- 
tion of ^■^l/T-^'-^l for all x shown by the solid grey curve. 
We show the procedure used to determine ^^xiocai in the 
case of a; = 0.07 shown as the solid black vertical and hor- 
izontal lines. In the case oi^^l/Tl^\ the horizontal black 
line determines what value oi^'^l/T^^l^ (solid grey curve) 
.(+) 



corresponds to ""^l/T^ , then the vertical line yields 
0.10 (±0.01) for the more metallic regions. A 



63 (+) 
local 

similar proceedure for ^'^l/T-^^~^ yields 
(±0.008) for the more insulating regions. The upper and 
lower frequency side therefore yield a consistent devia- 
tion ^^Axiocai = 0.028(±0.009) for x = 0.07 at 300 K 



We then determine 



63„( + ) 
■^local 

n(0) 



and 



63„(-) 
•^local 



for all X depend- 



ing on available ^ data and desired accuracy. We 

cannot determine ^^x\^^^i for x — 0.04 since we do not 
have ^^l/T^°A data below x = 0.04 at 300 K. The reason 
for this is that cannot be measured accurately for 

values > 7 (ms)~^. Furthermore, 300 K corresponds to 
the wipeout temperature T^qr for x < 0.035 |^, and 

we only determine ^^x\~^^i from our data for T > T^qr 
where full Cu signal intensity is observable. As for the 



limitations of we note that at 300 K, 

tends to merge to the same value of ^ 3 (ms)^^ for 



.(0) 



X > 0.09 |3^. This implies that one needs to measure 
1 /Ti """^ beyond the experimental uncertainties in order 

63^(+) 



local ■ 



to get a reliable estimate of a; 

The whole proceedure is repeated at different temper- 
atures, the final results of which are summarized in Fig. 
13 At different temperatures we find similar experimen- 



tal limitations in determining 



63 (+) 
local 



and 



63„(-) 
''local 



, how- 



ever, in cases such as a; = 0.07 where both can be deter- 



.(+) 

'local 



, 63 



Ax 



(-) 

local ' 



mined, we consistently find that ^^A; 
In the overdoped region x > 0.20 we find a compara- 
ble frequency dependence of ^^l/Ti across the NQR line 
to that of X = 0.16. However, as shown in Fig. |S1 we 

also find that ^^l/T^"^ starts to increase with increas- 
ing X for X > 0.20 at 300 K, and the same is true at all 
temperatures. As discussed in section II, we limit our 
determination of ^^xiocai to the region x < 0.16 and for 
clarity we separate the data in Fig. |Slfor x > 0.16 by 
using a dashed grey line. 

In Fig. ^1 we show the temperature dependence of 
^^1/Ti for X — 0.16 across the superconducting bound- 
ary Tc = 38 K (Fig. O at various positions across the 
NQR spectrum. Despite the fact that the spatial varia- 
tion ^^xiocai varies as much as 0.10 <^'^ xiocai < 0.22 at 
~ 100 K, all values of show a comparable frac- 

tional decrease below Tc where the superconducting gap 
opens. In the case of ^^1/T^~'> where '^^x'-'l^i = 0.10, 
should be around 25 K, however, there is already a large 
drop in 63i/r^(") by 30 K. In the case of 63i/y^(-i/io) 

where ^^x[~^l{^"^ = 0.08, Tc should be as low as 20 K, 

however shows a drop below 30 K. This 

shows that the superconducting transition is a genuine 
bulk phenomenon which effects all patches, regardless of 
the local hole concentration ^^xiocai- 



B. 



^UQ Results 



In Fig. Illf a') and (b) we show the temperature depen- 
dence of the resonance frequency at the CG of the NQR 
spectrum defined as (^^Vq) where k = (^4, B) for the A 
and B-lines, respectively. We define the first and second 
moments of the NQR spectrum as such 



^63 k 



ED.: 



63, ,fc / 63,,A; 



(12) 



where ^^v^ q corresponds to the jth data point of the 

observed NQR spectrum and ^Joq = \j2Ln2 = 1.177. 
For a Gaussian spectrum (found to be the case for x > 
0.07) the use of the prefactor y^oo reduces ^^Ai'q to the 
HWHM (half width at half maximum) exactly. 



8 



The (^^fq) data indicates two temperature regimes. 
Above Tq (black fit), we find a hnear decrease of I^^i^q) 
with decreasing temperature and a roughly constant 
value of {^^VqY In section IV we attribute the tem- 
perature dependence above Tq to thermal contraction of 
the lattice constants within the tetragonal phase, and in 
appendix A we also use this temperature region to de- 
termine the anti-shielding factors needed for the EFG 
simulation. Below Tq (grey fit), (^^i^q) increases with 
decreasing temperature for both A and B-lincs. The in- 
crease in this region is attributed to the increasing lo- 
cal orthorhombic distortions [s^ which in the case of 
X > 0.04 set in above the LRO structural phase tran- 
sition temperature Tst{< 515 K) according to neutron 
diffraction 0|. In the case of a; = 0.20, we find evidence 
for local orthorhombic distortion starting as high as Tq ~ 
350 K where diffraction results show no sign of LRO. In 
the section V we fit the data below Tq and deduce the 
degree of tilting of the CuOg octahedra away from the c- 
axis (defined as Oiocai), and we find that there is a sharp 
onset temperature for Oiocai only the case of a; = 0.0. 
For X > 0.0 we see a somewhat rounded transition into 
the local orthorhombic phase. We note that these local 
precursive effects to the LRO are consistent with pair dis- 
tribution function analysis of neutron powder diffraction 
data and XAFS (X-ray Absorption Fine Structure) 
analysis [i^ which find evidence for local lattice distor- 
tions in the temperature region above T^t, and also into 
the overdoped regime x > 0.20. 

The temperature dependence of the observed HWHM 
^^Avq of the NQR spectra, defined using Eq. ((T^ . is 
presented in Fig. 1121 We also reproduce the temperature 
regions above and below Tq, taken from Fig. ^2 as the 
black and grey curves respectively. In section IV we use 
the ^^Ai/q data for T > Tq to deduce a lower bound 
for the patch radius ^^Rpatch, and the corresponding an 
upper bound ^^Axpatchi'x I f"^ Rpatch) ■ 

IV. ANALYSIS 

A. EFG Background 

The goal of this section is to determine a static real- 
space model of the spatial variation in local hole con- 
centration ^^xiocai- We shall use a point charge lattice 
summation incorporating randomness effects to calculate 
the inhomogeneous distribution in the EFG which con- 
stitutes the observed ^^Cu NQR spectrum. To the best 
of our knowledge, such calculations incorporating ran- 
domness effects into the EFG calculation have not been 
reported. We shall deduce all parameters necessary for 
this task in a self-consistent way by using our pool of 
NQR data, and we shall successfully account for the en- 
tire ^'^Cu NQR spectrum, including the B-line, with one 
adjustable parameter i^^Rpatch) which defines the length 
scale of the spatial variation ^^xiocai- Once ^^Rpatch is 
optimized to fit the ^■^Cu NQR spectrum, we deduce an 



upper bound to the amplitude of the spatial variation 
across the plane ^'^ Axpatch{o^ 1/^^Rpatch), which we show 
is consistent with ^^Axiocai determined in a model inde- 
pendent way from section III. 

The distribution in the EFG is dominated by three 
mechanisms (a) the random substitution of donor ions, 
(b) the variation in local hole concentration ^'^xiocai 
over short length scales ^^Rpatch ^ 3.0 nm, and (c) 
the distribution in local lattice distortions. With- 
out a priori knowing the lattice distortions in (c), 
we attribute the whole NQR linebroadening to mech- 
anisms (a) and (b). This implies that the deduced 
^^Rpatch we present is a lower bound and correspondingly 

Axpatchifx. 1 /^'^ Rpatch) is an upper bound. However we 
shall demonstrate that ^^Axpatch Axiocah which im- 
plies that using mechanisms (a) and (b) alone is justified 
in the temperature range T > Tq. Below Tq where local 
orthorhombic distortions set in (Fig. Ill|l . we terminate 
our analysis of Axpatch- The reason for this is that we 
also see an increase in the NQR linewidth (Fig. I12fl below 
Tq in the case of a; < 0.07, which suggests that ignoring 
mechanism (c) above is no longer justified for T < Tq. 
In section V we make predictions about the magnitude 
of mechanism (c) by attributing the extra linebraodening 
for T < Tq to distributions in the lattice distortions. We 
note, however, that with the exception of section V, we 
ignore mechanism (c) in our analysis. 

There are two approaches for computing the EFG. 
The first is an ab initio approach psL liiL |45| which in- 
volves quantum chemistry calculations, and the second 
app roach makes use of the experimental data 0, 12M li^ 
I47ll48j to deduce all the necessary parameters in an em- 
pirical way. The ab initio approach has been extensively 
used to calculate the EFG at the ^^Cu site in La2Cu04 
psL 13 ^3 , and gives a consistent value for the observed 
resonance frequency (^^i^q) ~ 33 MHz jH to within 
uncertainties in the quadrupole moment ^'^Q. The ba- 
sic idea behind the ab initio approach is that since the 
EFG decreases rapidly ^ 1/r^ away from the origin of 
the calculation, the most significant contributions should 
be from local EFG components. One can then justifi- 
ably separate a cluster of ions in the immediate vicinity 
of the Cu nucleus whose contributions are calculated us- 
ing a full spin-polarized DF (density functional) or HE 
(Hatree-Fock) calculation, while the rest of the ions in 
the crystal are treated as point charges. The smallest 
realistic cluster is typically Cu06/Cu4Laio, where the 
central (CuOe)"^^ ionic cluster consists of 23 molecu- 
lar orbitals for each spin projection made up of linear 
combinations of 5 3d and 18 2p atomic orbitals. Such 
calculations are found to depend on the cluster size used, 
where the larger clusters are more reliable yet involve 
increasingly complex calculations. They also depend on 
whether one uses DF or HE, however both successfully 
predict that the molecular orbital with the highest en- 
ergy is the anti-bonding hybridization between the Cu 
3^3,2 _y2 orbital and the O 2pt, and 2py orbitals. The DF 
method, however, predicts [45j more covalent-like bond- 
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ing where the localized atomic spin density on the Cu is 
PCu = 0.67, while HF predicts more ionic-like bonding 
with pcu — 0.90. 

Before discussing our new approach for calculating 
the inhomogcneous distribution of the EFG, let us clar- 
ify exactly what inhomogeneous broadening implies in 
the context of ^'^Cu NQR. The linebroadening of the 
^^Cu NQR spectrum for x > 0.02 is dominated by the 
inhomogeneous distribution in the EFG. The inhomoge- 
neous broadening may be pictured as such: say there is 
a Cu nucleus j lying in a particular EFG with value Vj 
which resonates at frequency Vj, while a distant nucleus 
k sits in a distinct EFG environment Vk and resonates at 
I'k- As is predominantly the case across the sample, the 
separation between the two nuclei Y{j^k) is larger than 
the range ^af ~3 a of the indirect nuclear spin-spin cou- 
pling [43 where (,af is the correlation length of the anti- 
ferromagnetic fluctuations. In such cases, one may flip 
the j nucleus using an r.f. pulse and observe its echo 
without being effected by what the r.f. pulse does on 
the k nucleus. One then measures the ^^Cu NQR spec- 
trum by resonating at Vj and effectively counting the 
number of nuclei Nj in the sample with EFG values Vj, 
then changing the resonance frequency to Vk and count- 
ing the number of nuclei Nk, thereby building up the 
inhomogeneous NQR spectrum. 

As one approaches the undoped limit x < 0.02, the 
random effects become less and less significant, and the 
linebroadening becomes predominantly homogeneous in 
nature and dominated by the indirect nuclear spin-spin 
coupling [27j. As we show in appendix B, in the case of 
pure homogeneous broadening for x = 0.0, there is no fre- 
quency dependence in across the NQR spectrum, 
even in the orthorhombic phase T < Tgt- We note that 
the two distinct length scales £,af and ^^Rpatch which 
determine the linebroadening of the homogeneous and 
inhomogeneous NQR spectra, respectively, have qualita- 
tively different temperature dependences. As we shall 
show, ^^Rpatch mildly decreases with decreasing temper- 
ature, while ^AF increases with decreasing temperature 
[22]. Furthermore ^^Rpatch ^ £,af within the x and tem- 
perature range of interest. 

B. EFG Simulation 

In the experimental approach, the total EFG tensor 

'^tot'^^''^ (see Eq. (U))) along the orthogonal spatial di- 
rections (/3, 7) for the Cu site k = [A, B) is segregated 

into an onsite 2>d contribution Vg*-^''*^' and an offsite lattice 

contribution vj^^^'^^'^ as such 

T/(/3,7),fc _ /-fc T/(/3.'r) I /-fc T/(/3,7),fe /ION 

where we introduce the isotropic anti-shielding factor 
[27I Cm originating from the local electrons and the 
isotropic anti-shielding factor orginating from the lat- 
tice. We also leave the possibilty open that they depend 



on which Cu site k = {A, B) is being calculated. In ap- 
pendix A, we use a self-consistent approach to deduce 
that CLt = 27.7 and Ctd = 0.93, while Cu = 18-6 and 
— 0.82. We also show that the main principle value of 
the onsite EFG contribution V^j^ is negative while the lat- 
tice contribution VJ^'/j is positive, resulting in an overall 
negative value of V^^i . 

The anti-shielding factors must be determined exper- 
imentally, while in the ah initio approach they are, in 
a sense, already accounted for. The anti-shielding fac- 
tors are found not to vary significantly between differ- 
ent classes of high Tc cuprates Likewise, the anti- 
shielding factors deduced in a similar way for the planar 
oxygen site I'i, 46, J8J is also found to be consistent be- 
tween different classes of cuprates, which supports the 
use of the experimental approach. 

In accord with B. Bleaney and co-workers [Hj, the 
main principle value of the onsite contribution V^j^ arising 
from the 3dx2-y2 hole is taken to be 

Vi, = -^e{l-ia{r-,'), (14) 

where we have allowed for covalency in the form of f° 
which represents the fractional spin denisty on the four 
neighbouring planar oxygen 2po- orbitals. In order to 
be consistent with our results of ^''O hyperfine coupling 
analysis we use f° = 0.076 [T3 which favours more 
covalent-like bonding. Using ionic bonding with f° — 
results in an overall ~ 14 % increase in ^^Rpatch and 
corresponding ^ 12 % decrease in our theoretical esti- 
mate ^^Axpatch over covalent-like bonding. We note that 
f° = 0.076 deduced from the oxygen hyperfine couplings 
[13, 48] is consistent with Cu hyperfine coupling analysis 
using NMR ^3^1 where fcu (1 - 4/°) = 0.7 was found. 
It is also generally accepted that the ab initio DF calcu- 
lation which favors stronger covalency is a more realistic 
approach than the HF approach. 

The lattice contribution Vi^^^''^ to the total EFG is cal- 
culated using a point charge lattice summation with frac- 
tional point charges. We use the standard expression for 
the summation 

Viau^ = E 1^ (3^/3..^7,. - r'.Sp,,) lA. (15) 
3 

where the sum over j refers to the sum over the surround- 
ing point charges qj in the lattice a distance r-,- away. The 
charges are assigned as follows 

qcu = +(2-4/°)e 

go,, = -(2-2/°-x/2)e 

9o„p = -2e 

qha = +3e 

qsr = +2e, (16) 

where x is the nominal hole concentration. Op/ stands for 
the planar oxygen site, Oap stands for the apical oxygen 
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site and e(> 0) is the electronic charge (in emu). Note 
that we make use of the fuU difference in valency between 
the +^La and """^Sr ions. This implies that depending on 
the the geometrical distribution of the +^Sr ions in the 

lattice, vl|^^^'' takes on different values. In Fig. [T^ we 
illustrate an example of the random positioning of the 
^'^Sr ions surrounding the central Cu site from which the 
EFG is calculated. 

The standard approach in all earlier works, except for 
[T^ . has been to set q^a — Isr — — (3 — a;/2) which 
effectively bypasses any randomness effects and results 
in delta function spectra of the EFG. We start the EFG 
simulation by first taking ^^Rpatch = oo which uniformly 
places the donated holes x from the """^Sr ions with full 
mobility onto the planar O sites. The added holes reduce 
\qOpi \ in Eq. (|16|l . bearing in mind that there are two 
planar oxygens per unit cell. Placing the donated holes 
onto the planar oxygen sites as such is consistent with 
high energy spectroscopy studies js^l . 

The positions of the ions are allocated using an or- 
thorhombic cell (shown in Fig. I14|) which consists of 
4x (La2-£!;Sr^Cu04) formula units, each with the K2NiF4 
structure [Sj] . Using an orthorhombic cell will later per- 
mit us to incorporate the orthorhombic distortions in 
section V. The absolute values of the lattice constants 
[ao,bo,c] were calibrated at 295 K using our data from 
Fig. ^ The positions Yj{x,T) of the ions as a function 
of X and T were taken from a systematic X-ray powder 
diffraction study by P.G. Radaelh et al. The EFG 

calculation is naturally sensitive to inputs Vj (a:, T) , there- 
fore we use smooth interpolation of the values Vj{x,T) 
for different x in order to avoid unnessecary scattering 
in the output vlf^^\ We also use the same thermal ex- 
pansion coefficients 0111 eta = +1.45 • lO""'' K"^ and 
Qfc = -fl.42 • 10~^ K"-'^ for all x, consistent with cor- 
responding thermal coefficients found in La2-2;Ba2;Cu04 
p. 

The orthorhombic cell is then repeated in all space out 
to a radius of 50 A from the origin [0,0,0] of the EFG 
calculation. The procedure for the random placement of 
the """^Sr ions is carried out as follows: a different random 
number A with a flat probability distribution < A < 1 
is generated at each La site. If A > x/2 at a particular 
La site, a La ion with charge +3e is assigned to that site, 
while if A < a;/2, a Sr ion with charge -f2e is assigned. 
In equation form this gives the following prescription 

A > x/2 {-^) +3e 

A < x/2 +2e (17) 

at each La/Sr site. The factor 1/2 in the probability con- 
dition arises because there are 2 La/Sr sites per formula 
unit. Once this is carried out at each La site, we have 
achieved the particular random configuration k of """^Sr 
ions in the lattice, an illustration of which is given in 
Fig. 1131 In our formulation, we exclude the possibility 
of any +^Sr - "''^Sr clustering beyond probability theory, 
which is supported by the absence of any significant dif- 
fuse scattering in neutron diffraction experiments 54]. 



Once the random configuration k is determined, the 
summation in Eq. (|15|l is carried out and '^V^f^^'^ is di- 
agonalised and the principle value ''V^^tt stored into a 
vector Piatt- The calculation is computed in the tetrag- 
onal phase where the main principle axis of '^Vi^^^ ^ is 
found to lie primarily along the c-axis (i.e. ''VJ^^j is the 
main principle value). We note, however, that each con- 
figuration K, causes a slight rotation of the the main prin- 
ciple axis away from the c-axis. Strictly speaking, for 
each configuration k one should then go one step further 
and take the projection of the main principle value 
along the c-axis (i.e. along the average main principle 
axis). In the case of ^"^Cu NQR, however, the projec- 
tion only makes ~ 2 % difference in '^Vj^j, and can be 
neglected. 

Next, the whole lattice is re-randomized and a new 
random configuartion k' of +^Sr ions is determined, the 
lattice summation in Eq. H14|) is computed and the new 
value of V{g^^^ is stored into Piatt- This proceedure is 
re-iterated ^ 10^ times until the dimension of Piatt is 
sufficiently large to create a histogram spectrum of Piatt 
with ~ 50 bins across the lineshape, the results of which 
we show at 600 K in Fig. as the dashed grey lines. 
Note that we always calculate the EFG for the Cu site 
at the origin using a different random configuration of 
+^Sr ions for each run. An equivalent approach is to 
stick to one random configuration k for the whole sam- 
ple and calculate the EFG at each Cu site in the lattice 
thereby building up Piatt, however, this latter approach is 
computationally more intense. Since we are dealing with 
a purely random system, both methods are equivalent, 
therefore we use the former method. 

We notice three immediate features about the results 
in Fig. El The first is the overall decrease in CG (VJ^^j) 
with increasing x. This is a consequence of the change 
in lattice paramters and the decrease in IgOpJ with in- 
creasing X. The second feature is the change in HWHM 
which is found to increase with increasing a; as a result 
of the increased amount of quenched disorder in the lat- 
tice. The third feature in our lattice summation is the 
presence of the secondary peak known as the B-line. 

C. B-line 

The ^"^Cu NQR spectrum is known to show a secondary 
peak known as the B-line Generally, the origin of 
the B-line has been attributed to structural effects from 
the dopant ions P, , but the limited range of x in- 
vestigated in earlier studies did not allow unambiguous 
identification of the origin of the B-line. In fact, there 
has been a persistent claim that the origin of the B-line 
is more exotic in nature 3, 4, 43] and originates from an 
intrinsic response of the material to the presence of doped 
holes. Such a conjecture that the B-line is an electronic 
rather than a structural effect was based on the observed 
similarity between the A-B line splitting ^ 3 MHz for 
both La2-2;Srj;Cu04 and La2Cu04+5 systems USES- 



11 



It was thought that the same A-B hne sphtting between 
such distinct systems was clear evidence that the sec- 
ond site was a resuU of the presence of the doped holes 
themselves, independent of the means of doping. If this 
is the case, however, then it should equally well be ob- 
served in other 214 compounds such as La2-2;Baa;Cu04. 
However, a much larger A-B line splitting of ~ 6 MHz 
P, ^3 is observed in La2_a;Ba2:Cu04 where the +^Ba 
ion is known to cause much larger ionic size effect dis- 
order than the +^Sr ion, and which also stablizes 
an additional structural phase transition at low temper- 
atures |5^. We also recall that the ^■^Cu NQR spectrum 
in La2 .-„-i/Eu„Sr.T,Cu04 for ^^Cu isotope enriched sam- 
ples [lH shows a third structural peak (the C-line) whose 
fractional intensity is observed equal y. The C-line cor- 
responds to a Cu nucleus directly above or below a Eu^+ 
ion. 

In Fig. ^1 we present the experimental value of the 
relative intensity (f^^) of the B-line {Is) to the total 
{I A + Ib) defined as such 



rexp 
JB 



Ia 



(18) 



where both I a and Ib have been corrected for differences 
in spin-spin relaxation rates. Typically, ^^T'2,a/^^72^b ~ 
0.85 for both components of the spin-spin relaxation 
rate jZ^, similar to the ratio of ^^1/Ti. Our data in 
Fig. ^] shows that f'j^^ ~ x over a large range of 

0. 04 <x< 0.20, which is consistent with the EFG simu- 
lation where fs — x. According to the EFG simulation, 
the Cu B-sites (illustrated in Fig. 1141) correspond to Cu 
nuclei located at r = [riaa, rn,a, ricc] (where [rta, rtb, ric] are 
± integers ranging from zero to oo), which have distinct 
EFG values due to the presence of """^Sr ions located at 
positions r = [naa,nf,a, (ric ± 0.361)c]. In this interpreta- 
tion, the relative intensity of such Cu sites (/s) is then 
given by the concentration of +^Sr ions x, i.e. /b — x. 
This site assignment for the B-site Cu nuclei can natu- 
rally account for the observed relative intensity ~ x 

(Fig. USl). 

Other evidence that the B-line is primarily structural 
in origin is the similarity in the frequency dependence of 
63l/ri,s to ^^1/Ti^A (Fig. UJ. The similarity between in 
the frequency dependence of ^^l/Ti^k indicates that the 
length scale of the spatial variation in ^^xiocai is larger 
than the average +^Sr - +^Sr distance Isr = aj^/x, or 
equivalently the average B-site to B-site distance. Later 
we show that the deduced length scale ^^Rpatch ^ 3 nm 
is larger than Isr in the region of interest x > 0.02. Fur- 
thermore, as shown in Fig. |51 ^^xiocai across the B-line 
is consistent with ^^xiocai across the A-line. 

We also note that in Fig. H ^^Ti^a/^^Ti^b = 0.9-0.84, 

1. e. that ^^1/Ti^b is uniformly suppressed compared with 
^^l/Ti_A- Our systmetic study up to x = 0.20 showed 
that in the overdoped region x > 0.20, increases 
with increasing x. This implies that if there is hole lo- 
calisation in the vicinity of a Cu B-site as claimed in 
Ref. H S El, then ^^l/Ti^s at the B-site should be 



larger than ^^1/Ti^a in the case oi x = 0.16. Fig. 
El clearly shows that this is not the case, and further 
supports our interpretation that the B-line is primar- 
ily structural in origin. The most likely explanation for 
^^Ti,^/63Ti,b = 0.9 - 0.84 (and ^^T2,a/^^T2m) is that 
the hyperfine coupling constant |^'^^(Qaf)^| is 5-10 % 
lower at the B-site due to the increased lattice distortions 
in the vicinity of a "'"^Sr ion. It is known from XAFS [H^ 
that the +^Sr ion causes large structural distortions in 
its vicinity, which could well account for such changes in 
the hyperfine couplings. 



D. EFG 



We now proceed with our simulation taking the B-site 
assignment shown in Fig. E| We can separate the A 
and B-sites within the calculation by separating any ran- 
domly generated Vj^^^ value with a """^Sr ion in a position 
illustrated in Fig. 1141 and placing it into a new vector 
P^tt 1 while placing all other generated ''VJ^^j values into 
another vector P^^- This allows us to separate the A 
and B-lines as shown by the solid black lines in Fig. 1151 
and we attach the superscript k = (A, B) to '^Vi^tt- The 
first and second moments of the lattice contribution to 
the EFG are then computed as such 
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(19) 



In order to determine the total EFG defined in Eq. H13|l , 
we now determine the anti-shiclding factors for both the 
A and B-lines from the CG data as such 
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(20) 



where we have used Eq. H13I) and Eq. H19|l . along with 
the experimental data I^'^Vq) (see appendix A for more 
details). Eq. H20I) and the anti-shielding factors convert 
the CG of the EFG simulation to the CG of the experi- 
mentally observed data. 

We are now in a position to convert each EFG value 
("^/om) stored in Pl^^^ into a frequency as such 



2h 



(21) 



Using Eq. (|20|l and Eq. (|21|l . we can finally predict the 
intrinsic lattice lincbroadcning (Az^^^jJ as such 



2h 



Clatt ■ 



latt ■ 



(22) 



First, note that in Eq. and Eq. lfTO|l we have not at- 
tach a superscript k — {A, B) to AVJJ^jj. This is because 
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the EFG simulation predicts the same width AT^^^^ for 
both the A and B-Hnes. Second, note that we have as- 
sumed that the total width from the calculation Ai^/^^j 
is proportional to the intrinsic lattice width AVf^^^ with 
NO corresponding width from the onsite contribution, 
i.e. AV^^ = 0. We justify this by the fact that the ratio 
of widths between the A and B-lines is experimentally 
found to be 63Azy|/63A;y^ ~ 0.63 in the temperature re- 
gion T > Tq. According to appendix A, this ratio is 
closer to the ratio of the anti-shielding factors for the 
lattice contribution Cuitt/Ciau — ^■^'^ than for the onsite 
contribution (sd/^sd ~ 0.89. Note that the anti-shielding 
factors are deduced without any information about the 
widths, therefore taking AV^^ = is justified. 

Up to this point in the EFG simulation, we have taken 
care of the linebroadening due to the random placement 
of +^Sr ions in the lattice (mechanism (a)), however, 
we are still assuming that the holes are uniformily dis- 
tributed across the plane, i.e. ^^Rpatch = oo. Clearly this 
is not the whole picture since our calculation of Av^^^^ 
underestimates the observed linewidth ^^Avq by a fac- 
tor ~ 2, as can be seen in Fig. [21 



ion which satisfies the position 

r''^ ^ [(n„ + l/2)a, (nb + l/2)a,±0.139c] 



\Jn'l + nl < '^^Rpatch/a, 



(25) 



and in both cases, Ua and Ub are ± integers ranging from 
zero to oo. We do not a priori know which site di or ^2 are 
the donor sites, or if the donors come from a combination 
of c?i and c?2 • We therefore treat both or combinations of 
both as equal possibilities. Technically, the n.n. (nearest 
neighbour) c?2 is the closest hole donor to the central Cu 
nucleus when occupied by a +^Sr ion. The n.n. di site on 
the other hand has more bonding with the central Cu via 
the apical oxygens. The n.n di site has the largest effect 
on ^^VQ when occupied by a """^Sr ion and in particular 
gives rise to the B-site. 

Once we have counted all the N^^ satisfying either Eq. 
(|24|l or Eq. (|25|) for a particular random configuration 
of +^Sr ions k, we deduce the local hole concentration 
'^^^focai as such 



63 



^local = NUN' 



(26) 



E. 



^Rpatch Simulation 



The introduction of the spatial variation in local hole 
concentration ^^xiocai into the EFG calculation is an ex- 
tension of what we have described in arriving at the in- 
trinsic width Az^j^jj. We start by randomly positioning 
the """^Sr donor ions into the ^^La sites as described ear- 
lier. Next, we define a patch (i.e. a circle) of radius 
^^Rpatch around the Cu site at the origin, as shown in 
Fig. El We then determine the total number of Cu 
sites iV^p''*^'' within ^^Rpatch by multipying the area 
(T^^^Rlatch) patch covers by the areal Cu density 
(1/a^) as such 



^ ^patch 



(23) 



where we round N^p"*"*^ to the nearest integer. These 
]\lRpatoh gites within the radius ^^Rpatch will uni- 
formly share the donated holes from neighbouring +^Sr 
ions. Next, we determine how many donated holes Ng^ 
are to be shared within ^^Rpatch, where Ng^. itself is the 
random number to be calculated. As shown in Fig. 1141 
there are two distinct Sr/La sites (labeled di and 1^2), 
from which we can donate holes. This comes about be- 
cause there are two distinct Sr/La sites within the unit 
cell. For donation of holes from di, we accept a donor 
hole into ^^Rpatch provided there is a +^Sr ion which sat- 
isfies the position (in tetragonal notation) 



[uaCL, nta, ±0.361c] 



< 
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Rpatch/o, 



(24) 



relative to the origin. For donation of holes from d2, we 
accept a donor hole into ^^Rpatch provided there is a +^Sr 



We now carryout the EFG summation in Eq. (|14|l for the 
particular configuration of +^Sr ions k, where the planar 
oxygen charge qoj,, in Eq. 116|) is replaced by 



2/° 



63 



■^local 



/2) 



for all planar oxygens within the patch radius ^^i?- 
All other planar oxygen sites outside the radius ^^i?. 



(27) 

■patch • 
■patch ; 



including those on other CUO2 planes, are assigned the 
uniform hole concentration x. Using this scheme effec- 



tively ignores 3-d correlations of 



63 „ 



local 



from neighbour- 



ing Cu02 planes, however, such effects are small and of 
higher order. 

As before, we then convert '^V,'^'!l into a frequency us- 



latt 

into a vector P^. The 



ing Eq. (|21|l and we store 

position of the +^Sr ions are then re-randomized to ob- 
tain a new random configuration k', then a new value of 
i/'^ is deduced and stored into P^. This procedure is 
repeated ~ 10'* times until a sufficient number of data 
points exist in so as to create a histogram spectrum 
of the lineshape over ~50 bins. The CG of the spectra 
(Pj^) coincide with the experimental data {^^i'q} (Eq. 
(HOJ), while a new HWHM {Av^ ^^^J is determined as 
such 



63jyfc 



2. 



(28) 



We then repeat this whole procedure and calculate 
^^Rpatah ■various values of ^^Rpatck, the results of 
which are shown in Fig. El for the A-line. The results 
for the B-line are the same but with a uniform decrease 
^<<.../A<„... = CFatt/Ctti= 0.67) for all ^^Rpatch. 
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Also shown in Fig. El is how to deduce the best fit value 
of ^^Rpatch to the lineshape at 600 K. We match experi- 
mental values of Avq (shown by the grey horizontal line) 
to the calculation shown by the solid black curves, and 
extrapolate ^'^Rpatch shown by the grey vertical arrows. 

Note that according to the arrows, ^^Rpatch deduced 
using dopants from +^Sr ions at d2 is 15 % smaller 
than ^'^Rpl^tci^ deduced using dopants from +^Sr ions 
at di. Put in another way, when we do the calcula- 
tion for fixed ^^Rpatch, we deduce a ~ 15 % smaller 
value of Ai^*^ ^2 compared to A:/*^ . The reason for 

Ai^'^d2 < Aiy'^j^ can be explained by a correlation 

paten paten 

between the local hole concentration ^^xf^^^i on the pla- 
nar oxygens and the +^Sr donor ions in n.n. positions 
to the origin. The result of adding one hole from a n.n. 
+^Sr ion amounts to two separate shifts to the resonance 
frequency. The first shift is from the added hole itself 
which decreases of magnitude of the surrounding planar 
oxygen charge \qOpi\ and results in a large increase in 
the resonance frequency, regardless of the hole's origin. 
We shall define this shift as (5z/P'(> 0). The second effect 
of adding one hole from a n.n. """^Sr donor ion is from the 
n.n. +^Sr ion itself which results in a shift Sv^"^ , however, 
the sign of shift Si^^"^ depends on whether it is at ^2 or 

Let us imagine scenario (1) where the added hole 
specifically comes from a """^Sr ion in the n.n.n. (next 
nearest neighbour) di site. We do not use the n.n. 
di site in this case since this is the B-line which, ac- 
cording to the simulation, shows the same difference 
Av^a^ l^^v^d^ (^ 0.85) as the A-line. It can be 

^patch ^'patch 

shown that a +^Sr ion in the n.n.n. di site causes a small 
positive shift Sv'^^ > 0, therefore the total shift Sv^^^ in 
scenario (1) is given by (i5z/^'' -I- Sv'^'^) — <5:^*^^^(> 0) where 
both components are positive and enhance the total ef- 
fect. Next let us imagine scenario (2) where the """^Sr ion 
from which the added hole came from is specifically lo- 
cated at a n.n. d2 position to the origin. It can be shown 
that placing a """^Sr ion at the n.n. ^2 position causes a 
small negative shift (5i^'*^ < to the resonance frequency. 
The total shift Sv^^'^ in scenario (2) is therefore given by 
(^SiyP'' + Sv"^^) = <5i/*-^^(> 0) where the individual compo- 
nents counteract each other. Therefore, in scenario (1) 
the overall positive shift to the resonance frequency is 
enhanced while in scenario (2) it is reduced. After many 
random runs which include the independent random ef- 
fects from distant """^Sr ions, scenario (1) will lead to a 
larger overall spread, therefore Ai/'^rf^ < Av'^ 

'patch patch 



F. ^"^Rpatch Results 

In Fig. El we show the results for the best fit to the 
lineshape data at 600 K. Note that the observed asym- 
metry of the A and B-lines for x > 0.04 are well repro- 



duced within the model without any additional param- 
eters. The fits shown in Fig. El a-re deduced using di 
dopants with ^^Rpl^t^i^, however the same fits can be de- 
duced using d2 dopants with ^^Rp^i^i^ - The best fit values 

of ^^Katch and "^"^Rlltch are shown in Fig. ^&). In Fig. 
El we also show the data and fit for x — 0.02. As can 
be seen, however, the asymmetry in the experimentally 
observed line-profile for x = 0.02 is somewhat larger that 
our patch model can predict, hence we show the fit as a 
dashed line. It is interesting to note that in the range 
X < 0.02, 63i?^i^^^ and ''R^ltch tend towards the aver- 
age donor separation distance Isr = a/^/x, as shown in 
Fig. El If ^^Rpatch — hr, our modcl implies that each 
patch only covers one "^^Sr ion on average. In such a 
senario we no longer expect our model to account for the 
data. In fact, the critical region x = 0.02 which separates 
the spin- glas s phase and the anti-ferromagnetic phase for 
X < 0.02 (53 is known to show complex behaviour ^57j. 

We also note in Fig. El that the intensity of the B-line 
is slightly underestimated by the calculation compared 
with the experimental data. This is simply due to the 
fact that the spin-spin relaxation rate is experimentally ~ 
15 % smaller at the B-site together with the fact that we 
are measuring the data points using r = 12 /is. After cor- 
recting for differences in the spin-spin relaxation rate, the 
observed relative intensity of the B-line decreases slightly 
and is found to go as fj^^ ~ x (Fig. [TH|l . which is then 
consistent with the simulation. In order to fit the NQR 
spectrum in Fig. El we have attributed all of the extra 
line broadening to ^^Rpatch- In Fig. I19f b'l we compare 
the experimentally observed width ^^Avq with the cal- 
culated intrinsic width Av(^^^. 

In Fig. El we show the same plot as Fig. 1191 except at 
300 K, where the ^^Rpatch are deduced in a similar way 
to 600 K shown in Fig. El At 300 K, we see an effect 
not present at 600 K, namely that the observed linewidth 
^^Avq for X < 0.115 in Fig. I20r b') no longer decreases 
with decreasing x as it does at 600 K. Also shown in Fig. 
I20l is the onset of the LRO orthorhombic phase shown 
as the dashed lines through the data points. The result 
of the observed enhancement of the linewidths ^^Avq 
for X < 0.115 in Fig. I20f b) results in a large drop in the 
calculated ^^Rpatch shown in Fig. EJa). Comparing the 
X dependence in ^^Rpatch within the tetragonal phase at 
600 K and at 300 K for x > 0.115 seems to suggest that 
^^Rpatch is more or less x independent provided that the 
local lattice distortions do not dramatically enhance the 
observed NQR linebroadening. 

The effects of the orthorhombic distortions can also be 
seen in Fig. 1211 where we plot the temperature depen- 
dence of the average ^^Rpatch defined as such 

''RpatcH = {''R%tch Kltcn) /2. (29) 

We chose to plot the average value ^^Rpatch alone for 
clarity. This corresponds to a simulation where half of 
the donors come from c?i sites and half come from d2 
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sites. In Fig. [2] the black and grey lines correspond to 
regions above and below Tq taken from Fig. 1111 while 
the dashed black line shows the onset of LRO Tst j2^ . 
In the case of a; = 0.04 we clearly see a dip in ^^Rpatch 
below Tq ~ Tst which is most likely due to the fact that 
we have neglected distributions in the local orthorhombic 
distortions. In section V we make the reasonable assump- 
tion that ^^Rpatch for x = 0.04 continues to decrease lin- 
eaarly (defined as ^^R'patch) with decreasing temperature 
below Tq, similar to a; > 0.115 samples. Assuming that 
^^R'patch underlying length scale will allow us to es- 

timate the distribution in local orthorhombic distortions. 

We also remark that similar patch-by-patch models of 
disordered doping in the Cu02 plane have been indepen- 
dently used to account for the smeared density of states 
found in photoemission peaks as a result of a distribution 
in the superconducting gap ,58] . In such models the only 
free parameter used to account for the smoothness of the 
density of states is the length scale of the model which 
is determined by the in plane superconducting coherence 
length £^ab ^ 2.0 — 3.0 nm. It is interesting to note that in 
the case of x = 0.16, the lower bound ^^Rpatch ^ 2.5±0.5 
nm in the region of Tc is comparable to ^ab — 3.3 nm in 
thin films 59]. 

G. ^'■'Axpatch Results 

In the temperature region above Tq, we can also de- 
duce an upper boundary ^^Xpatch for the spatial variation 
^^xiocai without any further computation. We proceed as 
described above, but for each random lattice configura- 
tion K, we also store each value ^^xf^^^i (used to compute 
'^i^q) into a new vector Xj^^^^f . We then build up a vec- 
tor ^;oca/ randomly distributed local hole concentra- 
tions for each lattice configuration k. The mean of the 
set ^(oca; is just nominal value {Xf^^^j^ = x, while 
the effective HWHM ^^Axpatch (defined in an analogous 
manner to Eq. (|19|l 'l is found to be independent of which 
site k = {A, B) is computed. By deducing ^'^/S.Xpatch over 
a large range of the free parameters \x,^'^ Rpatch\ , we find 
our calculated values oi^^ /S.x patch fit well to the Binomial 
theorem 

CO xil — x) 

= 6315^ -^^(l^-l/^' (30) 

^•^Rpatch V 1^ 

where we have used the relation in Eq. H23|l for the sec- 
ond line. In Fig. IHlwc show the temperature dependence 
of ^^Xpatch using Eq. I© and Eq. (jH for ^^Rpatch de- 
fined in Eq. (EU. Using63i?di^^j63^d.^^j j^g^g^jgfgq^ 

(|29|l results in a small — (+)7 % variation in ^^Axpatch, 
respectively. The fact that ^^Xpatch xiocai suggests 
that (a) our patch-by-patch model is appropriate in de- 
scribing the spatial variation in ^^xiocah (b) our estimate 
of the intrinsic lattice broadening ^^Aviatt is correct in 



the temperature region above Tq, and (c) the static value 
^^Xpatch deduced from the lineshape analysis is consistent 
with the low frequency dynamic value ^^xiocai deduced 
using 63l/Ti. 



V. ORTHORHOMBIC LATTICE DISTORTIONS 
A. Local Order Parameter 

In the orthorhombic ph ase, the CuOe octahedra are 
tilted by an angle Oc < 5° |2J] away from the c-axis to- 
wards the &o-axis, and the primitive lattice constants in 
the plane are split [ho > ag)- We can calculate the ef- 
fect of these distortions to the EFG by placing the ions 
in an orthorhombic cell (see Fig . 1 1411 with appropri- 
ate orthorhombic modulations |2J, jsj and compute the 
changes to the ^^vq. However, we already know from 
Fig. ^2 that ^^VQ starts to increase above the structural 
phase transition temperature T^t deduced from neutron 
and X-ray diffraction results |2J] . These bulk probes de- 
tect LRO, whereas NQR is a local probe, therefore it is 
more appropriate to deduce the local orthorhombic dis- 
tortions from our ^^vq data, and then compare them with 
LRO results. 

The orthorhombic distortions cannot a priori be de- 
duced from the NQR data alone. We must first take cer- 
tain know relations from LRO results and apply them to 
NQR. In particular, we shall use the neutron diffraction 
results reported by P.G. Radaelli et al ^3] which demon- 
strate that certain measured orthorhombic parameters 
such as the orthorhombic splitting {ho ~ a-o), the octahe- 
dron tilting angle 9c{< 5°) which measures the departure 
angle of the Oap-Cu-Oap bond from the c-axis, and the 
"scissors" angle Oscis{< 90°) which measures the Op/- 
Cu-Op/ angle, are all related to powers of the underlying 
orthorhombic order parameter we define as 9lro{x, T). 

At low temperatures (< 70 K) where lattice fluctua- 
tions are small, mean-fleld theory give s a good account 
of the observed lattice distortions |2J]. In light of this, 
we fit the data in Ref. (2^ to a mean-field like form and 
deduced 

~ eLRo{x,T) 
esc^s ^ ■[BLRo{x,T)f 

{K~ao) ^ 6.1O-4.^^^.[0i«o(x,T)]^ (31) 

where all angles are in degrees. These set of equations 
provide a constraint for the orthorhombic distortions 
which determine the positions rj(x,T) of the jth ion in 
the lattice. We now assume that the LRO constraints 
in Eq. H31I) also apply to the local constraints which de- 
termine the local position Yj{x,T) of the jth ion in the 
lattice (i.e. eLR.o{x,T) 9iocaiix,T) in Eq. (EU)- 

The EFG lattice summation in Eq. H15|l is now carried 
out in a similar way as described in section IV but with 
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the orthorhombic distortions incorporated into the posi- 
tions of the ions. The positions rj {x, T) of the ions are de- 
duced using the absolute values of the lattice constants Co 
and c at 295 K from Fig. 0] along with their correspond- 
ing thermal coefficients aa and ac as before, however, we 
allow bo to vary according to the constraints in Eq. H31() 
and the unknown local structural parameter 9iocai{x,T). 
We position the La/Sr ions and the apical oxygen ions ac- 
cording to the tilting angle 6c = 6iocai{x, T) towards the 
bo direction, while the planar oxygens are positioned ac- 
cording to 9sc in Eq- and 6iocai{x,T). We now have 
all the ionic positions Vj (x, T) as a function of Oiocai {x,T). 
Just as in section IV, we then compute the EFG lattice 
summation and resonance frequency ^i'''^ (Eq. |2U) for 
the random configuration k and patch radius ^^Rpatch, 
however, the EFG lattice summation is taken using the 
tilted ionic positions specified by 9iocai{x,T). We store 
'^i/'^ into a vector Pg and re-randomize the lattice until 
Pg is sufficiently large to build a histogram spectrum. 
We define the new CG (Pg) incorporating finite 9iocai as 

The results of the simulation suggest that the increase 
in the CG due to 9iocai alone goes as 



(32) 



where c 



0.14 MHz/[deg.2] and c^(- 
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0.08 MIIz/[deg.^] are x and T independent within ~ ±15 
%. Note that NQR does not distinguish between ±6iocah 
but only the average square magnitude (9f^^^j^ across the 
whole sample. Another conclusion we deduce from the 
calculation is that after diagonalizing the EFG tensor, 
the main principle axis is found to be tilted an angle 
\9nqr\ away from the c-axis (where 9nqb^ — 0.9 • 9iocai), 
while the EFG asymmetry parameter rj remains small 
rj < 0.06. 

We are now in a position to deduce the average {9fgcai) 
from our experimental data {^^i^q) (Fig- 111(1 . In the 
spirit of the LRO results, we deduced the temperature 
dependence of {9f^^^i) (presented in Fig. I22II from the 



I^'^Vq) data assuming the paramteric form 



i^focal) 



Tq-T 
To 



2l3u 



(33) 



where {9fg^^i)'^ is the square averaged tilting angle at 
T = K, Piocai is the reduced temperature coefficient, 
and Tq is local onset temperature for the orhorhombic 
distortions. All three free parameters in Eq. H33|l depend 
on X and are optimized to fit {^^i^q) (in Fig. Ill|l . The 
best fit to the (^^q) are shown as the black and grey 
curves in Fig. 1111 and the corresponding temperature 
dependence of {9fc,cai) shown in Fig. In Fig. ESI we 

summarize the x dependence of best fit values {9fgcai)^ ^ 
Piocai and Tq . We compare these parameters with equiv- 
alent mean field parameters {9'^jiq) , /3 and T^t accord- 
ing to LRO |22,|2J|, shown as the solid grey lines in Fig. 



In the case of a; = 0.0, all three parameters deduced 
by NQR are consistent with equivalent LRO parameters 
deduced by neutron and X-ray diffraction. We find that 
Piocai — 0.34 according to NQR for x — 0.0, which is con- 
sistent with the LRO critical exponent f3 ~ 0.30 accord- 
ing to neutron and X-ray scattering for all x EOl - A 
similar critical exponent is also found in La2_2;Baa;Cu04 
where (3 = 0.33 ,53], which are all consistent with the 
theoretical predictions from a, 3d X ~ Y model. 

For X > 0.04 we find consistent values between 

i^focai) ^'^d {9°i^Ro) ■ III tlie case of Piocai and Tq for 
X > 0.04, however, we find qualitative differences with 
the LRO parameters (3 and Tst, respectively. The fact 
that Piocai ~ 1 for X > 0.04 indicates that locally, the 
transition into the orthorhombic phase is smeared out in 
temperature. The fact that Tq > 400 K indicates that 
the local onset temperature remains large and roughly x 
independent. 

In order to account for the differences between local 
and LRO structural results one must first recall that our 
NQR results measure the local square deviation 
while the diffraction results measure the coherent average 
{9lro) across the sample. One must also leave open the 
possibility that the orthorhombic structural transition for 
X > 0.04 is partly order-disorder and partly displacive in 
nature. In an order-disorder dominated structural tran- 
sition (as proposed by other local probes [4ll li^ in the 
region x > 0.15), 9iocai is argued to be finite and tempera- 
ture independent in the LRO tetragonal phase (T > Tst), 
while {9liio) — when averaged over large (> 10 nm) 
length scales. In such a scenario, one can picture cer- 
tain short length scale regions (<C 10 nm) which align 
along the positive bo direction where {9iocai) > 0, while 
an equal number of other regions align along the negative 
bo direction with {9iocai) < 0, therefore {9lro) = 0. The 
temperature dependence in an order-disorder type tran- 
sition does not involve a T dependence in the magnitude 
of |0;oca;| itself. Instead, the T dependence lies in the 
coherent orientation of the local regions towards a par- 
ticular direction, for arguments sake say along positive 
bo direction. The coherent average of the first moment 
{9lro) will then be finite and positive {{9lro) > 0), 
and will appear to increase with decreasing temperature 
as more and more local regions in the sample choose to 
align along the positive bo direction. 

In a displacive dominated structural transition (as ar- 
gued for X < 0.04), the magnitude of |0/ocaz| itself in- 
creases with decreasing temperature below Tst and all the 
local regions in the sample align along the same bo direc- 
tion. Our new NQR results combined with LRO results 
support a displacive structural transition for x < 0.04, 
while a mixture of displacive and order-disorder type 
transition occurs in the region x > 0.04. For x > 0.04 
below T < Tq, NQR indicates that locally the octahe- 
dra start to tilt in a displacive manner (i.e. the tilting 
angle of each octahedra increases with decreasing tem- 
perature below T < Tq and {9f^^^i) > 0). In the re- 
gion Tst < T < Tq, 9iocai continues to increase yet there 
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continue to be an equal number of regions with positive 
and negative tilting with respect to the direction, and 
therefore {9lro) = 0. In the region T < T^t, Oiocai con- 
tinues to increase with decreasing temperature, however, 
the local regions start to coherently align along one direc- 
tion (say along the positive bo direction), and (Olro) > 0. 
This interpretation is entirely consistent with other lo- 
cal probes 41, 42], but in addition we have deduced the 
onset temperature Tq for the local orthorhombic distor- 
tions. In the case of a; = 0.20 we find that Tq ^ 350 
K even though the LRO structural transition disappears 



B. Distribution in Local Parameter 

In all of our analysis of the ^'^Cu NQR linebroaden- 
ing thus far, we have ignored potential linebroadening 
from local lattice distortions. In the previous part of this 
section, we deduced the temperature dependence of the 
square average distortions {Ofocai) across the sample us- 
ing the CG data (^^i^q)- We shall now estimate the cor- 
responding distribution of Oiocai using the HWHM oi the 
NQR spectrum, ^^Ai/q. As discussed earlier, the three 
linebroadening mechanisms for the ^^Cu NQR lineshape 
are (a) the random placement of "^^Sr ions in the lattice 
(b) the patch-by-patch variation in hole concentration 
^^xiocai, and (c) the local lattice distortions. In section 
IV we used mechanisms (a) and (b) to account for the ex- 
perimental data, and deduced the lower bound ^^Rpatch 
for the spatial variation ^^xiocai- We shall now use all 
three mechanisms (a),(b) and (c) to account for the ob- 
served HWHM data ^^Av^ in the case oi x = 0.04. Note 
that we cannot a priori separate (a),(b) and (c), there- 
fore in order to obtain the contribution from (c), we must 
first make reasonable assumptions for (a) and (b). 

In Fig. 131 we plot the temperature dependence of 
the observed data ^"^Aj/q for x — 0.04 above and be- 
low Tq ~ Tst ^ 450 K taken from Fig. [El For T > Tq, 
we attribute all of the observed linebroadening to mech- 
anism (a) and (b) given by ^ ^ . For T < Tq , how- 
ever, we make the ansatz that the Av^, continues to 
increase only linearly, as shown in Fig. 1241 Ai/^, in 

Fig. EH results from using a different estimate ^^R'patch 
for the patch radius, shown as the grey dashed line in 
Fig. 12] for X = 0.04. The motivation for using '^^R'patch 
for X — 0.04 is that the temperature dependence fol- 
lows that of x > 0.115 where the overall amplitude of 

the orthorhombic distortions {{Of^caif ^1- (ESI) are 
smaller. 

In order to simulated the NQR spectrum, we now pro- 
ceed through the simulation as in section IV, however, 
we calculate the local hole concentration ^^Xi^^^i for a 
random configuration k using the patch radius ^^R'^atch 
(Fig. I21|) . instead of ^^Rpatch- Along with each simula- 
tion run K, we also independently choose a random value 



of |6'fgca/l taken from a trial PDF (probability distribu- 
tion function), examples of which arc shown in Fig. 1251 
The positions of all the ions in the random lattice con- 
figuration K are calculated given the value l^/'^o^a/l and 
Eq. (with Olro &iocai)- We then take the lattice 
summation in Eq. (|15|l using the the random positions of 
the ions with l^j'^j.^J, together with the random local hole 
concentration ^^Xi^^^^. The main principle value of the 
diagonalized resonance frequency '^v'' is then computed 
and stored into a vector Pg, then a new random lattice 
k' is re- generated along with a new value of |0;"j,aJ ac- 
cording to the trial PDF. The corresponding resonance 
frequency is calculated and stored into Pe . As in sec- 
tion IV, taking account of the random fluctuations of the 
main principle axis around the average NQR axis 9nqr 
(~ 0.9 • Oiocai) makes only minor ^ 2 % corrections to 
''z^*''. As before, we repeat this proceedure ~ 10^ times 
until Pe is sufficiently large to produce a histogram NQR 
spectrum with ~ 50 bins. Since ^^Rpatch (^^"^ ^ig. is 
already determined, the remaining free parameter in the 
simulation is now the PDF for |0;oca;|- 

The simulation is repeated using various trial PDF's of 
|6';oca;| until the the same simulation of the NQR spec- 
trum using ^^Rpatch from section IV is obtained. A good 
starting point for the PDF turns out to be a Gaussian 
distribution in the squared variable dfocai centered at the 
values (Of^cai) taken from Fig. ^ior x = 0.04. The 
width of the Gaussian distribution can then be optimized 
to best reproduce the observed NQR spectrum. After 
various trials, the final results of the optimized PDF's 
for |6'/oca/| (= \/ (^focai ) ^'^^ shown in Fig. |55l at various 
temperatures T < Tq. 



VI. CONCLUSIONS 

We have presented a detailed systematic study of the 
spatial variation in electronic states using ^^Cu NQR 
in ^"^Cu isotope enriched poly-crystalline samples of 
La2-2:Srj;Cu04 for 0.04 < a; < 0.16. By analysing the 
extent of the frequency dependence of ^^1/Ti across the 
inhomogeneous linebroadening of the ^^Cu NQR spec- 
trum, we determined the spatial variation in hole concen- 
tration ^^xiocai given by ^^xiocai = x AxiocaU where 
^^Axiocai was defined as the characteristic amplitude or 
extent of the spatial variaiton away from nominal hole 
concentration x. We showed that ^^xiocaii"^ x) is a ther- 
modynamic effect whose extent increases with decreasing 
temperature below 500-600 K and reaches values as large 
as ^'^Axiocai/x ~ 0.5 in the temperature region > 150 
K. Furthermore, we showed that the spatial variation 
^'^xiocai is an intrinsic effect in La2-a;Srj;Cu04 by con- 
firming the same extent of the inhomogeneity in our poly- 
crystalline samples to that in high-quality single crystals 
hj. To the best of our knowledge these reports, together 
with those in Ref.'s [l^ IT^ . are the first of its kind to 
show the temperature dependence of the intrinsic inho- 
mogeneity in La2_2;Sr3;Cu04, or any other system with 
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quenched disorder. In Ref. the results of the inho- 
mogeneous electronic state ^"^xiocai determined using ^^O 
NMR in high-quality La2-a:Sra:Cu04 crystals are shown 
to be consistent with the ^^Cu NQR result ^^xiocai- We 
believe that this is a powerful demonstration of the va- 
lidity of our analysis in this work, as well as in 0, ■ 
Moreover, consistency between ^"^O NMR and ^'^Cu NQR 
results indicate that the spatial variation in the spin and 
charge channels are correlated across most of the Bril- 
louin zone. 

We have argued that there are two essential ingredi- 
ents for the spatial variation ^^xiocaiij^ x). They are 
as follows: (a) that the +^Sr donor ions are positioned 
with random probability throughout the lattice and (b) 
that there is a short length scale for the spatial varia- 
tion ^^xiocai of the order ^^Rpatch ^ 3.0 — 4.0 nm. Us- 
ing (a) and (b) above we succesfuUy fit the inhomoge- 
neous ^'^Cu NQR spectrum (including the B-line) using 
a patch-by-patch model for the spatial variation ^^xiocai 
with the patch radius ^^Rpatch as the only adjustable 
parameter. We assumed that no extra ^^Cu NQR line- 
broadening from the lattice distortions existed, therefore 
the calculation for ^^Rpatch resulted in a lower bound. As 
shown in Ref. |l3j |. however, our lower bound estimate 
^^Rpatch is shown to be consistent with ^"^ Rpatch deduced 
from the inhomogeneous ^^O NMR spectrum, which jus- 
tifies our assumptions used to determine ^^Rpatch in the 
temperature region T > Tq. Within the EFG sim- 
ulation, we also deduced an upper boundary to the 
spatial variation in electronic states ^^Xpatch given by 
^^Xpatch — X zL^^ Axpatch, where the extent of the spa- 
tial variation ^'^Axpatch was shown to depend on ^^Rpatch 
as such ''^Axpatch oc if'^Rpatch (i-e. '^^Rpatch = oo im- 
plied the homogenous doping limit ^^x patch = x). Within 
the temperature region T > Tq we found the consistent 
result ^'^Xpatch xiocai, which justified our EFG simu- 
lation and short length scale model of the inhomogeneous 
electronic state. 

^^Rpatch was deduced assuming no clustering of the 
+^Sr ions in the lattice beyond probability theory, how- 
ever, we note that even in the case of clustering over 
length scales 3-4 nm, the presence of a short length 
scale is still needed to account for the observed spa- 
tial variation in electronic states. We also argued that 
the short length scale implied that any experimental 
probe which averages over spatial extents larger than 
^^Rpatch ^ 3.0 — 4.0 nm results in an underestimate of 
the extent of the inhomogeneous electronic state. NQR, 
however, is a local probe, which makes it ideal for de- 
tecting short ~nm variations in the electronic state. Our 
findings have a strong impact on the general, but now 
incorrect, view that doped holes are homogeneously dis- 
tributed in the Cu02 planes of La2-xSra;Cu04. 

We have also compared our local orthorhombic struc- 
tural distortions results to LRO structural results taken 
from neutron and X-ray diffraction [2^ 0| techniques 
which probe coherence phenomena over length scales 
larger than ten's of nm '25'|. In the case of x = 0.0 we 



found good agreement between NQR and LRO results, 
however, for x > 0.04 we presented evidence that the lo- 
cal orthorhombic distortions started at elevated temper- 
atures Tq{> 400 K) above the onset of LRO Tst- We also 
found that the local distortions were somewhat smeared 
below Tq compared with x = 0.0. Our local structural re- 
sults were found to be consistent with a mixture of order- 
disorder and displacive transitions in the region x > 0.04, 
consistent with other local probes [111 |42| . Our NQR 
studies at elevated temperatures have now quantified the 
onset of the local orthorhombic distortions at Tq{> 400 
K) in the region 0.04 <x< 0.20. 

It is known that in the a: < 1/8 region of the phase di- 
agram, the temperature dependence of the wipeout frac- 
tion ^^F{T) (defined as the fracton of unobservable Cu 
nuclei that are wiped-out) is somewhat tailed jS. iTll 
below Tnqr in both La2-xSr2:Cu04 and its Eu and Nd 
co-doped materials. In the case of La2-xSra;Cu04 for 
X — 0.04, the onset temperature for Cu wipeout is rather 
high at Tmqr — 300 K, yet half the Cu nuclei are still 
observable (i.e. ^'^F{T) ~ 1/2) at ~ Tmqr/2, implying 
that T^QR does not signify the onset of a collective phe- 
nomenon for X < 1/8. This may indicate the formation 
of a small number of local moments in the Cu02 plane 
which cause the wipeout of the NQR signal for neigh- 
bouring Cu sites, as in a classic spin-glass. 

Rather than invoking a spin-glass type model, however, 
we can now naturally account for the tailed wipeout for 
X < 1/8 using the temperature dependence of the grow- 
ing spatial variation in electronic states ^^xiocai- In the 
temperature region T^qr ^ 300 K for x = 0.04, the char- 
acteristic local hole concentration ^^xiocai for insulating 
patches is inferred to be as low as ^^xiocai ^ 0.015 (Fig. 
0. Recall that ^H/Ti for x = 0.02 at - 300 K is already 
too large for the Cu signal to be observable, therefore the 
more insulating patches for x = 0.04 at 300 K begin to 
wipe-out. As one cools the x — 0.04 sample below 300 
K, there are a growing number of patches with lower and 
lower hole concentration which are also wiped-out, lead- 
ing to an gradual increase in the wipeout fraction ^^F{T) 
with decreasing temperature. The fact that ^'^xiocai only 
gradually fans out with decreasing temperature naturally 
explains the tailed nature of the wipeout fraction ^^F{T) 
below Tj^QR in the underdoped region x < 1/8. 

In Eu and Nd co-doped La2-a;Sra;Cu04 for a; < 1/8, 
the same tailed feature in ^^F(T) is found together with 
an inflection point T^"(^J|'=(< Tnqr) below which ^^F{T) 
increases rapidly with decreasing temperature. In the 
case of Eu and Nd co-doping with x < 1/8, T^''q^'^ 
coincides with the onset temperature T^tt of an addi- 
tional structural phase transition with tetragonal sym- 
metry known as the LTT (low temperature tetragonal) 
phase and the onset of charge ordering Tcharge- In the 

case of Lai.6-xNdo.4Sra;Cu04, Tcharge — Tltt ^ 70 K 

[g^ . This coincidence suggests that the glassy slow- 
ing of the stripe modulation below T^q^^(~ Tcharge) 
is accelerated by the tetragonal symmetry compatible 
with the charge stripe symmetry |15|. ^'^Cu isotope en- 
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riched samples of Eu and Nd co-doped La2-a:Sra:Cu04 
are found to have comparable frequency dependence of 
^^l/Ti across its inhomogeneous ^'^Cu NQR spectrum. 
In the Eu and Nd co-doped cases, therefore, we can nat- 
urally account for the tailed wipeout in the temperature 
region T^q^'^ < T < Tnqr by the growing inhomo- 
geneity in electronic states with decreasing temperature, 
while for T < Tj^q^^(~ Tcharge), the wipeout is dramat- 
ically accelerated due to the pinning of charge stripes by 
the lattice symmetry. Although there is no LRO into 
a LTT phase in La2-a;Sra;Cu04 which can stablize the 
charge stripes, there is evidence for an incipient struc- 
tural transition into the LTT phase , which may sug- 
gest a similar phenomenon is occuring in La2-xSra;Cu04. 



APPENDIX A: ANTI-SHIELDING FACTORS 



In this appendix we deduce the anti-shielding factors 
|27l | C;^(( and CidJ where k = {A,B). We first assume 
that both anti-shielding factors are isotropic and inde- 
pendent of X and T, but are allowed to vary between the 
A and B-lines. We produce a cross-plot (shown in Fig. 
HBJa)) of the resonance frequency at the CG of the A-line 
(^^i^q) [x,T) (Fig. inia)) against the CG of the lattice 



contribution to the EFG ( V, 
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{x,T) (Fig. 113). Like- 
wise, we produce the same cross-plot for the B-line shown 
Fig. I26r b). Note that we only use data in temperature 
region above Tg (Fig. Ill|l where local orthorhombic dis- 
tortions are absent. 

Then we deduce the linear fit using the standard em- 
pirical form [IEIStI: 



calculation to experiment are the A's in Eq. IjAip . Like- 
wise, when accounting for the linebroadening in section 
IV to deduce ^^Axpatch, it is also the A's which ulti- 
mately matter, not the C's, therefore our end result for 
^^Axpatch and ^^Rpatch are independent of what value of 
^^Q we use. 



APPENDIX B: MAGNETIZATION RECOVERY 

In this appendix we quantify the observed multi- 
exponential recovery of the magnetization M{t) in the 
case of a: = 0.16. The fact that the recovery is multi- 
exponential implies that is distributed at a fixed 
frequency on the NQR line. We can naturally account 
for the distribution in ^^1/Ti at fixed frequency by re- 
calling that the ^"^Cu NQR line has an underlying intrin- 
sic lattice broadening (Aj/j^jJ which will tend to smear 
out the frequency dependence of ^^1/Ti across the NQR 
line. Smearing out the frequency dependence of 
results in a distribution of ^^1/Ti values at a fixed fre- 
quency. Using the the multi-exponential recovery data, 
we estimate the extent of the underlying lattice broaden- 
ing Avf^l^ = 0.62(±0.07) MHz in the case of x = 0.16 at 



the A-line. We find that AiyifJ^^ is consistent with the lat- 
tice broadening Ai^^^^ — 0.49 MHz found independently 
using our point charge calculation (section IV), which 
also adds weight to our EFG simulation of the random 
lattice. 

The experimentally observed nuclear recovery M(t) 
can be fit to Eq. (HIJ, where M(0), M(oo) and 
are free parameters of the fit. The bare recovery is then 
defined as I{t) where 
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where k — (A, B) and all quantities are overall positive. 
The best fit to the data shown by the solid black lines in 
Fig. indicates that A^^ = 75.3 MHz and A^„ = 16.0 

64.6 



MHz/[emuxlO~^''] for the A-line, and that Af^ 
MHz and A;^„ = 10.7 MHz/[emux 10""] for the B-line. 

Using Eq. H20(l we express the A's in terms of the 
unknown anti-shielding factors as such 
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Using /° = 0.076 ^ ,^^Q = -0.16 barns jp, and the 
bare value (r^^) = 7.5 a.u. for free cupric ions, we 
find that for the A-hne = 27.7 and = 0-93, while 
for the B-linc (iatt = 18-6 and Cm = 0.82. These set of 



results are consistent with previous estimates [23, |43 . 

Converting the above analysis using ^^Q = —0.211 
barns 13(1 instead of ^^Q = —0.16 barns 31] results in a 
uniform 24.2 % decrease in our estimates of the C's. How- 
ever, note that what really matters in converting from 
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In the case of single exponential recovery, I(t) is then 
found to be a straight line on a semi-log plot, an example 
of which is shown in Fig. I27l where we present I{t) in the 
case X = 0.0 at 475 K for various positions across its NQR 
spectrum. As discussed earlier, the x = 0.0 spectrum is 
homogeneously broadened, and therefore shows no fre- 
quency dependence of ®^1/Ti across its line (i.e. I{t) has 
the same slope at different frequencies), nor does it show 
any sign of having a distribution at a fixed frequency. 
We also note that the data in Fig. |^is taken at 475 K 
which is the orthorhombic phase, therefore we can also 
rule out the possibility that the orthorhombic distortions 
are the direct cause of observed frequency dependence in 

63l/Ti. 

In Fig. |2Slwe show the bare recovery I{t) at similar po- 
sitions across the A-line for a; = 0.16 at 100 K. First note 
that the slopes in the recovery get steeper towards lower 
frequency. This corresponds to the frequency dependence 
in 63i/Ti where ^^^l/T^^' >63 i/yW >63 i/t[+}. Sec- 
ond, note that I{t) shows finite curvature on a semi-log 
plot. This implies that there is a distribution in recover- 
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ies I(t) which can be fit to the general form 
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Fitting the recoveries to a single exponential form (i.e. 
taking aj =1 and l/Tij = 1/Ti in Eq. ||B2l)) results 

in the force fit values ^^1/T^% and ^^1/T^^2- 

These force fit values represent the average of the under- 
lying distribution in ^^l/Ti ji. 

The dashed lines in Fig. |2H1 show the best fit result 
from a sum of exponentials l/Tij whose distribution co- 
efficients Gj are presented as the dashed lines in Fig. 1^ 
The distribution coefficients aj in Fig. I29l are determined 
as such 



(B3) 
(B4) 




where Eq. l|B3(l is used for 1/Tij < Pti and Eq. l|B4|) 
is used for 1/Tij > Pti- The three free parameters used 
to fit I{t) are therefore the peak of the distribution Pti 



and the two widths Ap and A'^'. When optimizing the 
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three parameters in Eq. ljB4|) . we also use the constraint 
that the CG of the distribution (shown as the solid ver- 
tical lines in Fig. I29() coincide with the force-fit value 
^^l/Ti.A (presented in Fig. deduced with a single 
exponential fit. 

We then deduce the second moment of the distribu- 
tions in Fig. in an analogous manner to Eq. H12|) 
and Eq. (HHJl, and we find that the HWHM oi the distri- 
butions are Al/T^^2 = 0-33 (ms)-\ Al/T^'"] = 0.27 
(ms)~i and Al/r/+) = 0.20 (ms)-i. Note that the 
HWHM oi the distributions are ~ 15 % of average value 
which implies that accurate recovery data I{t) 
is needed down to I{t) ~ 0.001. 



We can now get an estimate of the underlying intrin- 
sic lattice broadening Ai/fjj.^ by converting the HWHM 
Al/Ti^A deduced above into equivalent frequency widths. 
In order to do so we take a local derivative /3,y of the un- 
derlying frequency dependence of at a frequency 
v defined as such 
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where is the interpolated value of ^^1/Ti a shown 

in Fig. 1201 as the solid line. Using ^^1/Ti,a in units of 
(ms)^^ and v in units of MHz, we find that /S*^^) — 0.62, 
/3(o) = 0.39, and = 0.29, from which we deduce the 
HWHM (in MHz) 
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0.53 



A1/t/"]//3(") =0.69 



Aj.(+|^ = A1/tW//3^+^ = 0.68 ^ (B6) 
for the lower half intensity, center and upper half inten- 
sity of the spectrum, respectively. Finally, we deduce 
the average value AvfJ^^ of the three estimates in Eq. 
(|B6|) to be Avf^^^ = 0.62(±0.07) MHz. This is consistent 
with Au{l^^ — 0.49 MHz deduced independently using 
the point charge calculation, as shown in Fig. 12111 

We note that our whole analysis in this appendix is 
derived using fixed t — 12 /is conditions. As discussed in 
section II and HI, however, there exists a small r depen- 
dence on the underlying distribution in the spin-lattice 
relaxation rate. Such effects could account for the 25 % 
difference between the lattice width AvfJ^^ derived from 
Eq. (|B6p and the point charge calculation Avf^^^ in sec- 
tion IV. 
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FIG. 1: (a) Frequency dependence of ^^I/Ti.a and "^^I/Ti.b 
across the A and B-hnes of the ®^Cu NQR spectrum of ^'^Cu 
isotope enriched La2-a;Sri,Cu04 for x — 0.20 (x), x = 0.16 
(A), X = 0.115 (•), X = 0.07 (o) and x = 0.04 (a). All data 
is taken at 300 K. Solid curves in (a) are a guide for the eye. 
(b) ^"^Cu NQR spectrum of the A and B-lines where the same 
symbols as part (a) are used. Solid curves in (b) are calculated 
fits using patch-by-patch model of spatial variation in local 
hole concentration ^^xiocai- The only free parameter used 
in the fits is the lower bound for the patch radius ^^Rpatch = 
2.1 — 3.1 nm (using di dopants, see later) in order of increasing 
X, respectively. 



FIG. 4: X dependence of the lattice constants at 295 K de- 
duced by X-ray diffraction along the Uo (o), bo (x), and c (•) 
axes using orthorhombic notation. Lines are a guide for the 
eye. 



FIG. 5: Temperature dependence of ^^l/Tl"i at the CG of the 
A-line for x = 0.20 (solid grey line), s = o!l6 (A), a; = 0.115 
(•), x = 0.09 (♦), X = 0.07 (o), X = 0.04 (a), x = 0.02 
(o), and X = 0.0 {+). Sohd black curves are a guide for 
the eye. All data (in this Fig. and in Fig. |SJ| are taken 
above the ^^Cu wipeout temperature Tnqr using a fixed pulse 
separation time of r = 12 ^s. 



FIG. 2: Frequency dependence of ^^1/Ti,a and ^^1/Ti,b (•) 
across ^^Cu NQR lineshape (o) for x = 0.115 at 300 K taken 
from Fig. ID We define ^^i/T^C-i/io) ^ ''''l/T^, "'l/Ti'"!, 

and at the lower one-tenth, the lower half, the CG 

(center of gravity), and upper half intensity position of the A- 
line, along with ®^ 1 /T^'^g at the CG of the B-line. Curves show 
fit using patch-by-patch model with a lower bound for the 
patch radius of ^^Rpatch = 1-6 nm (dotted curve), ^^Rpatch ~ 
2.6 nm (solid curve), and ^^Rpatch ~ oo (dashed curve, whose 
HWHM is given by Ai/^t)- 



FIG. 6: Temperature dependence of ®^l/r^^"^ at the CG 
of the B-line where the same data symbols from Fig. 13 
are used, together with the black solid lines which repre- 
sent (^e-^^ 1/Ti°i) data taken directly from Fig. with a 
uniform scaling factor e, where e = [0.87, 0.87, 0.90, 0.84] for 
X = [0.20,0.16,0.115,0.07], respectively. Grey dashed line 
represents ^'^l/T^'^l for x = 0.20. 



FIG. 7: Temperature dependence of ^^I/tI^^^^^ (+), 

""^l/T^ (.), '^''1/T^% (black curve alone), and (°) 
with nominal hole concentration x given in each panel. Grey 
curves are '^^l/T^^ for x = 0.00, 0,02, 0.04, 0.07 ,0.09, 0.115, 

and 0.16 where ^'^l/Tj'^^ monotonically decreases with in- 
creasing X. 



FIG. 8: X dependence of ^'^l/Tl'J (•), "^^l/Tli (x), and 

(°) at 295 K where the vertical dashed lines connect 
the data for each x. Grey solid curve shows interpolation of 

for ^1 while dashed grey line shows overdoped 
regime. The solid black horizontal and vertical lines illustrate 
how to extract ^^xiocai iu the case of a; = 0.07 according 
to the upper frequency side of the ®^Cu NQR 

spectrum defined as ^^ajj^caii ^'^'i according to ^''^1/t[ 2 
Shirane, J. Phys. Soc. Jpn. 69, 851 (2000). ^he lower frequency side defined as '''a:,'-'^,. 



FIG. 3: Temperature dependence of field-cooled Meissner sig- 
nal in 15 Oe for poly-crystalline La2-a^Sra;Cu04 with x — 0.20 
(x), X = 0.16 (A), X = 0.115 (•), X = 0.09 (♦), x = 0.07 (o), 
and X = 0.04 (a). Lines are a guide for the eye. 



FIG. 9: Temperature dependence of the local hole concentra- 
tion ^^xiocai (grey filled •) deduced from the upper frequency 
side of the A-line using '^^l/T^'^2 C^^loLi > and also from 
the lower frequency side using ^^l/T^'2 C^x\~^^i < x), where 
X is indicated in each section and shown as the grey horizon- 
tal lines. Solid black lines show reflections of (grey filled •) 
data through x lines. Also shown is ^''^xiocai deduced using 
^^1/T^^2 and ^^1/T^~J at the B-line (■), together with the 
upper boundary ^^Xpatch (x) deduced from calculated fit to 
the '''^Cu NQR spectrum using a patch-by-patch model for 
the spatial variation. Representative error bars for ^^xiocai 
are also shown. 
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FIG. 10: Temperature dependence of ''^l/Ti for x = 0.16 
across the superconducting boundary at Tc = 38 K (shown as 
grey verticle line) at various positions across the ^"^Cu NQR 
line including ^^/t'^^/"') (+), ''H/t[;2 (•), ''1/^^ (x), 
'^^1/t['^} (o), and ^^l/xfl (▼). Black curves are a guide for 
the eye. 



FIG. 16: X dependence of the fractional B-hne intensity /s 
defined as /s = Ia+Ib ' ^^^'^^ ^'^'^ have been corrected 
for differences between the spin-spin relaxation rates at the A 
and B-sites, respectively, at 300 K (+) and 600 K (x). Solid 
line shows prediction from model /s = and an illustration 
of the B-site is specifically shown in Fig. 1141 



FIG. 11: Temperature dependence of CG of the ^'^Cu NQR 
spectrum (a) (^^Vq^ at the A-line, and (b) {^^i^q) at the B- 
line for x = 0.20 (x), x = 0.16 (A), x = 0.115 (•), x = 0.07 
(o), X = 0.04 (A), X = 0.02 (o), and x = 0.0 (+). Black 
curves show fit above Tq where ^^Xpatch is deduced, while grey 
curves show fit below Tq where local orthorhombic distortions 
are deduced. Dashed black line shows onset of orthorhombic 
phase according to LRO jS^j. 



FIG. 13: Section of the Cu02 plane illustrating an example of 
the random placement of ^'^Sr ions. The corner of each square 
represents a Cu site and the central Cu site (o) is taken as the 
origin of the EFG calculation. The (grey •) symbols repre- 
sent the planar O sites. The ^^Sr donor ions (•) immdediatley 
above or below the Cu sites are shown, corresponding to the 
donors di (see later). Also shown is a typical patch with a 
radius of ^^Rpatch ~ 3.0 nm (=8 a) from the origin. For this 
particular random configuration k, the local hole concentra- 
tion within the patch is '^xiocai ~ 0.05. 



FIG. 17: Calculated value of HWHM Au^^^^^^ (MHz) for 

the A-line as a function patch radius ^^Rpatch (nm) for dif- 
ferent X shown in each part using the "'"^Sr dopant ions di 

^^flpatch.di (°)' ^^'^ dopants A;y^ ^2 (•)• Solid black 
curves through data are guides for the eye. Grey solid hor- 
izontal line shows experimental data ^"^Auq at 600 K while 
dashed black line shows intrinsic lattice width Au^ reached 
in the limit ^^Rpatch ~ 00. 



FIG. 14: Illustration of the A-site Cu nucleus (•) and B- 
site Cu nucleus (grey filled o) determined according to ad- 
jacent positions of "'"^La ion (grey •) and ^'^Sr ion (grey o), 
respectively. Also shown are methods of doping into the cen- 
tral plane using di ^^Sr sites (black curves) or d2 '^^Sr sites 
(grey curves). An orthorhombic cell [ao,bo,c] is used, while 
dashed grey line illustrates Cu-Cu lattice spacing a in tetrag- 
onal phase. 



FIG. 15: Results of the EFG simulation incorporating the 
random positioning of """^Sr ions in the lattice. Grey dashed 
lines indicate spectrum of the principle value Vf^^t of the EFG 
in [emu xlO~^*] at 600 K for various x shown in each part. 
Black lines show the same data but with the A and B-sites 
separated (V,^j^, where k = {A,B)) within the calculation. 

Short grey vertical lines indicate CG position (v^^^^ used to 
calculate anti-shielding factors. 



FIG. 12: Temperature dependence of HWHM of the ^'^Cu 
NQR spectrum (a) Auq at the A-line, and (b) ^''^ Auq at 
the B-line for x = 0.20 (x), x = 0.16 (A), x = 0.115 (•), 
X = 0.07 (o), X = 0.04 (a), and x = 0.02 (o). Black and 
grey curves indicate temperature region above and below Tq 
respectively. 



FIG. 18: ^^Cu NQR lineshape at 600 K for x = 0.16 (A), 
X = 0.115 (•), X = 0.07 (o), X = 0.04 (A), x = 0.02 (o) and 
X = 0.0 (grey region). Solid lines are best fits to the spectra 
using the patch-by-patch model and di dopants. The lower 
bound to the patch radius Rpatch ~ 3.0 — 4.0 nm used for 
the solid line fits are shown in Fig. 1191 
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FIG. 19: (a) Lower bound estimate of the patch radius 
^'"^Rpatch (in units of nm on left axis and in multiples of lat- 
tice spacing a on right axis) used to fit spectra in Fig. 1181 at 
600 K using the ^^Sr donor ions di (o) or d2 (•) as donors, 
(o) corresponds to "'"'^Sr-^^Sr separation Isr = a/y/x. (b) Ex- 
perimentally deduced HWHM ^'^Auq (■) of spectra shown 
in Fig. 1181 along with calculated intrinsic lattice broadening 
^'^latt foi" fhs A-line. All lines are a guide for the eye. Also 
shown is Avf^l^. (x) for x = 0.16 deduced in appendix B. 

FIG. 20: Same plot and symbols as Fig. EH but at 300 K. 
Dashed black lines corresponds to orthorhombic structural 
phase for x < 0.115 according to LRO [23l . All other lines 
are a guide for the eye. 



FIG. 21: Temperature dependence of lower bound estimate 
to the patch radius ^^Rpatch (defined according to average in 
Eq. deduced from the width of the ''^Cu NQR spectrum 

shown in Fig. [TH for s = 0.16 (A), x = 0.115 (•), x = 
0.07 (o), and x = 0.04 (a). Lines are a guide for the eye. 
Grey lines indicate the temperature region T < Tq where 
local orthorhombic distortions are present (Fig. IIH . Dashed 
black line corresponds to LRO structural temperature Tst l23 . 
Dashed grey line for x — 0.04 corresponds to ''^R'patch (see 
text) for T <Tq. 



FIG. 22: Temperature dependence for the square of the local 
structural parameter (or equivalently, the square of the CuOs 
octahedron tilting angle) (Of^^^i^ (in deg.^) for x = 0.20 (x), 
X = 0.16 (A), X = 0.115 (•), X = 0.07 (o), x = 0.04 (a) and 
X — 0.0 (-f ) deduced from the data in Fig. 1111 



FIG. 23: X dependence for the coefficients of the local struc- 
tural parameter (Ji'iocai) including (a) maximum square value 

{^locai)" (deg.^), (b) reduced temperature coefficient Piocai, 
and (c) onset temperature Tq (K) for x — 0.0 (o) and x > 0.0 
(•). Dashed lines are guides for the eye. Solid grey lines are 
LRO patameters (a) {el^of (deg.^) [H, (b) /3 [EEi^] and 
(c) Tst 22] . 



FIG. 24: Temperature dependence of experimental HWHM 
^^Av^ for = 0.04 (a) taken from Fig. [TH where black 
and grey lines show temperature region above and below 
Tq = Tst = 450 K. Also shown is Ai/^, ^ ^ (o) for T < Tq de- 
termined using R'patch taken from Fig. 1211 The remaining 
linebroadening due to orthorhombic distortions alone, Au^ 
(x), are also estimated (see Ref. [6l|l. 
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FIG. 25: Probability Distribution Function for the magnitude 
of local tilting angle |&iocai| for x = 0.04 at various tempera- 
tures T < Tq. At each temperature, the PDF's shown result 
in the best fit to the NQR lineshape assuming ^'^R'patch from 



from Fig. [TTIa) against the calculated EFG at CG (Vi^^^j 
taken from Fig. llSl with x and T as implicit parameters, both 
for the A line, (b) Same cross-plot as in (a) but for the B-line. 



FIG. 27: Time dependence of bare recovery I{t) deduced 
from magnetisation M{t) (see Eq. lB2t for x — 0.0 at 475 
K (< Tst), taken at various positions across NQR spectrum 



including '>H/Tl"i (grey x), ^'l/n+J (o) and ""'l/n^J (•). 



Solid line is best fit to a single exponential recovery. 



FIG. 28: Same plot and symbols as Fig. HTjbut for i = 0.16 
at 100 K. Solid lines are force fits using a single exponen- 
tial recovery, while dashed lines are best fits using multiple 
exponential recoveries whose distribution is shown in Fig. 1291 



FIG. 29: Distribution coefficients aj of l/Tij (dahsed lines) 
deduced from best fit to recoveries I{t) shown in Fig. 1281 for 
X = 0.16 at 100 K. Solid vertical lines show both force fit 



values '^^1/T|7 ^^l/Tl"^ (grey line) and ^^l/T^^^ in order of 



increasing frequency shown by arrow which coincide with CG 
of distributions. 



FIG. 30: ''^Cu NQR spectrum of x = 0.16 at 100 K (•) along 
with force fit values ^'^l/Ti (o) at various positions across the 
line. Grey lineshape shows fit using patch-by-patch model 
with a lower bound for the patch radius of ^^Rpatch = 2.7±0.2 
nm, while grey dashed lines shows intrinsic linebroadening 
(with HWHM of A/zf^tt) according to EFG simulation (sec- 
tion IV) . Black lineshape shows intrinsic linebroadening (with 
HWHMoi l^vTatt) deduced from recoveries in Fig. 1281 at the 
A line. Solid curves through (o) are interpolations between 
data points. 
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